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INTRODUCTION 


In the early 1970's, the Langley Research Center of The Na- 
tional Aeronautics and Space Administration initiated a combined 
experimental and analytical investigation to study the aerodyna- 
mic characteristics of irregular planform wings (also sometimes 
referred to as cranked wings (ref. 1) or double- delta wings 
(ref. 2). For this study, the planforms were referred to as wing- 
fillet combinations with the inboard more highly swept portion of 
the planform being defined as a fillet. 

The early phase of the study was directed toward improving 
the aerodynamics of the space shuttle orbiter (e.g. , ref. 3), al- 
though the general long-range goals are applicable toward improved 
design of aircraft as well as certain advanced aerospace vehicles. 
The benefits to be derived from the use of fillets with selected 
planforms include linearization of the subsonic lift- curve slope 
to high angles of attack. With regard to the space shuttle orbi- 
ter design, the improved lift at the angle of attack specified for 
landing allowed for either reduced landing speed or reduced wing 
planform area for specified mission return weight. In addition, 
proper tailoring of the wing-fillet combination allows lineariza- 
tion of the curve of pitching moment against angle of attack to 
angles for high lift; thus, trim penalties on both lift and per- 
formance are reduced. Although these subsonic benefits might be 
favorable, the question arose as to what effect a near-optimum de- 
sign would have on the desired hypersonic trim angle and stability 
requirements (dictated by cross-range or heating constraints). 
Since both subsonic and hypersonic conditions were the two prime 
areas of concern in the application of wing-fillet combinations, 
the overall study was designated the Subsonic-Hypersonic Irregu- 
lar Planform Study (SHIPS) . 

With regard to the overall SHIPS program, the objectives of 
the study are to generate an experimental data base from low sub- 
sonic to hypersonic speeds accounting for secondary effects of 
Reynolds number, airfoil section, leading-edge radius and sweep 
as well as planform geometry; to provide an aerodynamic predic- 
tion technique for irregular planform wings based on these exten- 
sive wind-tunnel results; and to provide empirically determined 
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boundaries to serve as design guides regarding linearized lift, 
pitch, and realistic longitudinal center-of-pressure locations 
as functions of Mach number. 

The present report is an element of the overall SHIPS pro- 
gram and presents the results of an investigation to develop a 
prediction technique for the low-speed static aerodynamic charac 
teristics in pitch of a class of low- aspect ratio irregular plan 
form wings for application in preliminary design studies of ad- 
vanced aerospace vehicles. 

The presentation is organized in the following manner. 

•The scope of the investigation is discussed first 
to provide a basis for understanding the goals and 
the technical approaches used. 

•The experimental data base is described in enough 
detail to support certain judgments that were made 
during the course of the investigation. 

•A presentation of the many geometric parameters 
used in the investigation and the equations used 
to generate them is provided. 

•The results of an evaluation of previously exist- 
ing prediction methods are presented and discussed. 

•Efforts made to develop additional correlations of 
test data to help formulate new prediction methods 
are described and results presented. 

•The development of the basic set of prediction 
methods for lift, drag, and pitching moment of 
irregular planforms having NACA 0008 airfoils is 
described. 

• The analyses accomplished to account for the secon- 
dary effects of Reynolds number, airfoil thickness, 
airfoil thickness distribution and airfoil camber 
are illustrated and modifications made to the predic- 
tion methods are presented. 
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•Comparisons between predicted lift, drag and pitch 
ing moments and test data are presented. 

•Concluding remarks are presented. 


•The effects of Reynolds number on the lift-correla 
tion parameter are described in an Appendix. 
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A BASIC = *W 
TRUE 


a. c . 


a. c 
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7oC 


b 


c 



Aspect ratio 

Aspect ratio of basic wing 
Aspect ratio of irregular planform wing 
Aspect ratio of inboard wing panel 
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Aerodynamic- center location, in percent 
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Drag coefficient 
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LIST OF SYMBOLS AND NOMENCLATURE (Cont’d.) 


'D 
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'D 
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J D 
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'D 
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'D . 
min 


'D. 


MISC 


'f . 
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TURB 


Base-drag coefficient 
Cross-flow drag coefficient 
Friction-drag coefficient 
Form-drag coefficient 
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Drag due to lift 

i 
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coefficient 
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C^ - Lift coefficient produced by Spencer's empiri- 

SPENCER cal method for irregular planforms 

Lift coefficient produced by WINSTAN empirical 
method for sharp leading edge double delta 
planforms 

Lift- curve slope evaluated near zero lift, 
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Design lift coefficient for two-dimensional 
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Cjjj - Pitching-moment coefficient 
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LIST OF SYMBOLS AND NOMENCLATURE (Cont'd.) 


Cm. 25c = “ Pitching -moment referenced to quarter chord 

CM-.25CBAR location of mean geometric chord 


Cma 


Slope of pitching- moment variation with angle 
of attack evaluated near zero angle of attack 


W* 0,1, 2, 3, 

4, 5, 6, 7 


Slopes of linear segment in pitching-moment 
prediction method 


c =MGC= CBAR - Mean geometric chord, meters (ft) 


c e ff = e WF ” Mean geometric chord of irregular planform, 

meters (ft) 


£ ref 


Value of mean geometric chord used as reference 
length for pitching-moment coefficient 


F 


1 ’ 


Functions in White-Christoph skin-friction 
analysis 


FF 

FR 


Form factor in minimum-drag prediction 

Body fineness ratio for arbitrary crossection 
body 

Fuselage station location of aerodynamic 
center, inches 


f 

f 
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Function in skin- friction analysis 

Stall progression factor used in lift predic- 
tion 


h b 


Height of rectangular part of body cross section, 
meters (ft) 
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LIST OF SYMBOLS AND NOMENCLATURE (Cont’d.) 


H f 


Control surface hinge factor used in minimum- 
drag analysis and prediction 

IF 

- 

Interference factor used in minimum-drag 
analysis 

K 

- 

Admissible surface roughness (equivalent sand- 
grain roughness) diameter, cm (inches) 

K 1 

- 

Function in admissible roughness prediction 

L 

- 

Characteristic length for calculation of Rey- 
nolds number in skin-friction analysis, m (ft) 

l B 

- 

Overall body length, meters (ft) 

LER? 

c w 

- 

Leading-edge radius at mean geometric chord 
of basic-wing planform, meters (ft) 

LER r 

- 

Leading-edge radius of tip chord, meters (ft) 

L N 

- 

Forebody length, meters (ft) 

l ref 

- 

Reference length for pitching-moment coefficient 

l 

- 

Overall length of irregular planform, meters (ft) 

M 

- 

Mach number 

MGC = MAC = 
c = CBAR 

- 

Mean geometric (aerodynamic) chord length, 
meters (ft) 

R B 

- 

Radius of semicircular portion of body cross- 
section, meters (ft) 
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LIST OF SYMBOLS AND NOMENCLATURE (Cont'd.) 



Reynolds number based on leading edge radius 
and velocity evaluated normal to leading edge 




Lifting surface interference factor used in 
minimum-drag prediction ; 


R, 


N 


Reynolds number 


R N L 


Reynolds number based on characteristic length 
of aircraft component 


r N/L 


Unit Reynolds number, 1/meter (1/ft) 


R NLERc w 


Reynolds number based on leading-edge radius 
of mean geometric chord of basic wing 


r nler t 


Reynolds number based on leading-edge radius at 
tip 


Rm * 
w mm 

R Nmax 


minimum Reynolds number 
Maximum Reynolds number 
Standard Reynolds number 


defined in text 
for use in pre- 
diction method 



Wing-body- interference factor used in minimum- 
drag calculation 


"R" 


Leading- edge- suction ratio 



Plateau value of leading-edge- suction ratio 
from correlation of basic wing data 



Plateau value of leading-edge- suction ratio for 
irregular planforms 


tlTJ M 

R (a) 


Leading-edge- suction ratio as function of angle 
of attack 


9 



LIST OF SYMBOLS AND NOMENCLATURE (Cont'd.) 



'R 


«*) 


1 , 2 , 3 , 
4 , 5 , 6, 7 


Value of leading-edge-suction ratio in various 
regions of linear segmented representation 


r 


Leading-edge radius of airfoil, meters (ft) 


r/c 


Leading-edge radius as decimal fraction of 
airfoil chord 


r 


S BASE 


Recovery factor used in skin friction analy- 
sis only 

Base area contributing to base drag in mini- 
mum drag prediction, Sq. meters (Sq.ft) 


REF 


TRUE 


’W 


WET 


Reference area for aerodynamic coefficients 
(prediction met 
meters (Sq. ft)) 


(prediction methods use Spgp = S^^p, Sq, 


Total planform area of irregular planforms , 
Sq meters (Sq ft) 

- Wing area of basic wings, Sq. meters (Sq.ft) 

- Wetted area of body, Sq. meters (Sq. ft) 


B 




s 




Wetted area of inboard wing panels, Sq. meters 
(Sq. ft) 

- Wetted area of outboard wing panels, Sq meters 
(Sq- ft) 

Wing semispan, meters (ft) 

Slenderness ratio = s/l 
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LIST OF SYMBOLS AND NOMENCLATURE (Cont'd. ) 


Taw 

T 


t 

t 


t/c 

W 

W ! 


w 


B 


X 


X 


c.g 


X. 

1 


X 


(x/c) 


a . c . 


Adiabatic wall temperature 
Free-stream temperature 
Function in skin-friction equation 
Airfoil thickness - meters (ft) 

Airfoil thickness as fraction of chord 
Basic-wing planform 

Inboard wing panel of irregular planform 
Outboard wing panel of irregular planform 
Body width, meters (ft) 

Longitudinal distance from apex of planform, 
meters (ft) 

X location of pitching-moment reference point, 
meters (ft) 

X distance from apex to break in leading-edge 
sweep 

Longitudinal distance from body nose or wing 
leading edge to point of transition to tur- 
bulent flow, meters (ft) 

Chordwise location of aerodynamic center as 
decimal fraction of wing chord 
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LIST OF SYMBOLS AND NOMENCLATURE (Cont’d.) 


(x/c) 


t/c 


max 


Chordwise location of airfoil maximum thick- 
ness as fraction of chord 


x 


y 

a 

tt BREAK 


X distance from apex of planform or panel to 
leading edge of mean geometric chord 

Spanwise distance from root chord of planform 
or wing panel to mean geometric chord 

Angle of attack, degrees or radians 

Angle of attack for upper limit of primary- 
slope region of pitching moment, degrees 


a. 


a STALL 


Angle of attack for initial flow separation 
used in Chappel’s wing flow separation analysis, 
degrees 

Limit of linear lift in WINSTAN cranked wing 
analysis , degrees 


“ 1 , 2 , 3 , 4 , 
5 , 6,7 


Angle-of-attack boundaries in suction ratio 
and pitching-moment predictions , degrees 


a. 


‘'ref * ref 

4ref> ^ref 


$ 

8 a 


( 1 ) , 2 , 3 , 
4,5 


Angle-of-attack boundaries for reference Reynolds 
number condition in suction-ratio prediction, de- 
grees 

Prandtl Glauert factor: p = 

Incremental value of angle-of-attack boundary 
due to a change in Reynolds number, degrees 


'LvOR 


SLE 


Incremental value of lift coefficient caused 
by vortex flow from sharp wing leading edge 
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LIST OF SYMBOLS AND NOMENCLATURE (Cont'd.) 


Aa.e . break “ Change in aerodynamic center location at angle 

of attack for upper limit of primary slope re- 
gion of pitching moment 



Increment of lift coefficient above linear 
theory 



Increment of lift coefficient for round leading 
edge 


Ac l S le 


Increment of lift coefficient for sharp leading 
edge 


Ax 

A(sa ( )CosA w j . 



Aa (3 _2) 

%4-3) 

A «(5-4) 


Fictitious distance ahead of transition in mixed 
laminar -turbulent flow prediction for skin fric- 
tion, meters (ft) 

Incremental value of change in angle -of -attack 
boundary function due to a change in Reynolds 
number, degrees 


Incremental value of lift correlation parameter 
(note three different types of increments are 
discussed in the text) 

Incremental value of angle of attack between 
boundaries Ot^ and > degrees 

Incremental value of angle of attack between 
boundaries OC ^ and Oi^ » degrees 

Incremental value of angle of attack between 
boundaries OL 5 and OL^> degrees 
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LIST OF SYMBOLS AND NOMENCLATURE (Cont'd.) 


Vv = ETAB 

D 

®LE 

A 



A TE 


Nondimens ional spanwise location of break in 
leading-edge sweep 

Compliment of leading-edge-sweep angle , degree 
Sweep angle, degrees 

Leading-edge -sweep angle of fillet, degrees 

Leading-edge -sweep angle of basic wing and 
outboard wing panel, degrees 

Trailing-edge- sweep angle, degrees 


a c/2 

A c/4 

A t/ C 
u t c max 

Q 


Sweep angle of half chord line, degrees 
Sweep angle of quarter chord line, degrees 

Sweep angle of maximum thickness line, degrees 
Taper ratio of basic wing planform 

Correlation function for suction ratio in 
plateau region 



Correlation function for effect of Reynolds 
number on angle-of -attack boundaries defining 
regions of linear variation of suction ratio 
with angle of attack 

Longitudinal stability derivative evaluated 
near zero angle of attack 
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LIST OF SYMBOLS AND NOMENCLATURE (Concluded) 


) Slope of suction -ratio curve in each region 

( 2 ), 3 , 4 , 

5,6,7 


CODES 

Configuration code used in many figures 

a F Aw .Airfoil 
Example 80 45 .0008 

Computer Procedure Code 

RlT - Air Force version of Aeromodule Computer Proce- 
dure (Ref . 8) 

X6I - SHIPS Aerodynamic Prediction Procedure for Ir- 
regular Planform Wings 


Note: Quantities are presented in the International System of 

Units (U.S. customary units in parenthesis). The work 
was performed using U.S. customary units. 
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SCOPE OF THE INVESTIGATION 


The scope of this investigation as originally conceived was 
based on the analysis of wind-tunnel test data from a single para- 
metric experimental investigation of 35 planforms having NACA 
0008 airfoils plus a few configurations having NACA 0012 airfoils 
or sharp leading-edge double-wedge airfoils. The nominal unit 
Reynolds number range of the test data was from 6.56 million per 
meter (2 million per foot) to 26.25 million/meter (8 million per 
foot). The test data base consisted of 131 pitch runs. 

Additional tests were tentatively scheduled to provide data 
to higher Reynolds numbers on a limited series of planforms in- 
cluding variations of airfoil sections. Data from the additional 
tests were to be considered when and if they became available. 

The basic objectives of the investigation were: 

(1) To evolve empirical methods from the SHIPS 
experimental data base for the prediction 

of first and second order subsonic lift, drag, 
and pitching- moment characteristics of ir- 
regular planform wings of moderate to high 
thickness ratios having application on pos- 
sible advanced aerospace vehicles. 

(2) To provide correlating parameters and simple 
predesign charts as well as a rapid and ef- 
ficient computer program for quick evalua- 
tion of new configurational concepts. 

The proposed technical approach to meet these objectives had 
four basic elements. 

(1) Existing prediction methods would be examined 
first to evaluate their applicability. 

(2) Where needed, correlations of the SHIPS ex- 
perimental data would be accomplished to 
formulate improved methods. 
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(3) The basic prediction methods would be de- 
veloped from the data available from the 
initial test at the maximum unit Reynolds 
number for the configurations having 
NACA 0008 airfoils. 

(4) Modifications to the basic methods to ac- 
count for Reynolds number effects and air- 
foil section effects would be sought de- 
pending on the scope of data available 
during the course of the investigation. 

In fact, five additional tests were accomplished by NASA 
which increased the test data base to 452 pitch runs as de- 
scribed in the next section. The additional tests were essen- 
tial to meet the objectives of the investigation. 

Prediction methods were sought for the following elements 
of the static aerodynamic characteristics in pitch. 

• Lift- curve slope near zero lift ^C. 

•Nonlinear lift increment ' 

•Low angle-of -attack stability derivative 

•Aerodynamic center location - (x/c) a c 

• Angle of attack for stall 

•Pitch- up /pitch- down boundaries 

•Drag due to lift - C D 

L 

•Variation of leading-edge-suction ratio with angle 
of attack - M R (a) " 

It was assumed from previous experience that existing 
methods of predicting the minimum-drag coefficient were suf- 
ficiently accurate to meet the objectives of the investigation 
if properly applied to the irregular planforms. 

The magnitude of the analysis task is illustrated by the 
fact that more than 4000 curves were plotted during the course 
of the investigation. Small programmable calculators and desk 
top computers were beneficial in manipulating and plotting the 
large mass of data. 
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EXPERIMENTAL DATA BASE 


The scope of the experimental data base is described in this 
section. The planform study (ref. 4) consisted of 35 planforms 
illustrated in Figure 1 in which the geometry is shown normalized 
with respect to the root chord of the basic wings. The detailed 
geometric characteristics of these planforms are presented in 
Tables 1 and 2 for the initial series of models. In essence, the 
planform families started with five basic tapered planforms having 
leading-edge sweeps of 25, 35, 45, 53 and 60 degrees. The irregu- 
lar planforms were generated by adding various fillets of increased 
leading edge sweeps up to a maximum of 80 degrees. The wing area 
(S^) and aspect ratio (A BASIC . } .265) of the basic planforms were 
constant, and the spanwise location of the intersection of the fil- 
let leading edge and the basic wing leading edges was constant 
(T)g =0.41157), The airfoil chordwise thickness distribution was 
constant across the span. 

The models were constructed such that each planform was a 
separate model which was attached under a minimum cross section 
balance housing as shown in Figure 2. Nose fairings and constant 
cross section aft extensions were fitted to the balance housing to 
produce a "minimum body" for each wing- fillet combination. 

Data were supplied from six tests as described below. 

(1) Test ARC 086-12-1. Ames Research Center 12-foot pres- 
sure tunnel. Planform Study - small models. 

Complete planform matrix (35 planforms) with 

NACA 0008 airfoils 

3 planforms with NACA 0012 airfoils 

A w = 25°; a f = 25°, 60°, 80° 

4 planforms with 8 percent thick modified double 

wedge airfoils 

A w = 35°, A p = 35°, 80° 

A w = 60°, A f = 60° , 80° 
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Boundary-layer transition strips on wings and body 
nose for most runs. 

Mach number = 0.30 

Nominal unit Reynolds numbers = 6.56, 13.12, 19.69, 
26.25 million per meter (2, 4, 6, 8 million per foot) 

NOTE: some planforms were tested at only the two 

or three highest unit Reynolds numbers. 

a Range from -2° to 26° 

SrF-F = total planform area - STRUE 
LREF = MGC of total planform =c^p 
X CG = 0.25 MGC 

(2) Test 8-ft TPT-780 Langley Research Center Eight-Foot 

Transonic Pressure Tunnel. Inboard Panel Alone Tests - 
small models. Zero suction wing test. 


AIRFOIL A f = 65°, 
NACA 0012 
NACA 0008 x 

THIN MODIFIED DOUBLE 
UNI T REYNOLDS NUMBER 

MACH NUMBER 


70°, 75°, 80° 

(x) 

XXX 

WEDGE A „ = 45°, A_ = 70° 

9.84 (11.48) (16.4) x 10°/meter 

3.0 (3.5) (5.0) x lCT/foot 

0.3 (.6) (.8) 


NOTE: Conditions in ( ) run only for NACA 0012 

model with A_ = 80° 

F 

Boundary- layer transition strips on body nose and 
wing surfaces 

^Range: -2° to 21.5° 


^REF = ^ ota l planform area of irregular planform* 

Lj^, f = Mean geometric chord of irregular planform* 

Xcg = 0.25 mean geometric chord of irregular planform 

*For fillet-alone tests, planform area was that of 
Wing II plus fillets. 
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(3) Test LTPT-255. Langley Research Center Low-Turbulence 
Pressure Tunnel. Airfoil Thickness Study to High Rey- 
nolds Numbers . 


Limited Matrix of Planforms - small models 


AIRFOIL SECTION A = 45°; A = 
NACA 0008 
NACA 0012 
NACA 0015 

THIN MODIFIED DOUBLE WEDGE 
UNIT REYNOLDS 

NUMBERS 13.12, 19.67, 26.25, 
4, 6, 8, 

MACH NUMBER .30 .30 .30 


45 ( 

65°, 

o 

o 

r^. 

, 75 

0 

1 , 

80 

0 

1 


X 


X 

X 


X 



X 

X 

X 

X 


X 



X 


X 

X 


X 



X 



X 


X 



32, 

.81, 39 

.37, 

45. 

93 

X 

10 6 /M 


10, 

12, 


14 

X 

io 6 /ft 


.28 

.26 
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Boundary -layer transition strips on wings and body nose 
for most runs. Complete range of unit Reynolds numbers 
for each configuration. 

a Range from -2° to 26° 

S REF = basic win S area " S W 

Lp^,p = Root chord of irregular planforms -C^ 

= 70 percent of root chord 

(4) Test LTPT-262. Langley Research Center Low -Turbulence 
Pressure Tunnel. Additional airfoil section studies to 
high Reynolds numbers. 


Limited Planform Matrix - small models 


AIRFOIL SECTION A„ = 45°; 

A p = 65°, 

o 

o 

75°, 

o 

o 

00 

NACA 0008 

X 




NACA 0012 

X 




NACA 0015 

X 




NACA 6 5 A0 12 


X 

X 

X 

NACA 65,412 


X 

X 

X 
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UNIT REYNOLDS 

NUMBERS 13. 12 , 19.67, 26.25, 32.81, 

4 6 8 10 

MACH NUMBER .30 .30 .30 .28 

Free transition for all runs. 

a Range from -2° to 36°+ 

S KEF = % 

Lref = Root chord of irregular planforms -C^ 

X CG = P er cent of Root chord -C^ 

(5) Test LTPT-266. Langley Research Center Low-Turbulence 
Pressure Tunnel. More additional Airfoil Section 
Studies to High Reynolds Numbers . 


39.37, 45.93 x 10 6 /M 
12 14 x 10 6 /Ft 

.24 .21 


Limited Planform Matrix - small models 


AIRFOIL SECTION A T7 

W 

NACA 0008 
NACA 0012 
NACA 0015 
NACA 0004 
UNIT REYNOLDS 

NUMBERS 13.12, 19 
4 

MACH NUMBERS . 2 
(.3) 


45°; A p = 60°, 65°, 

X 

X 

X 

X XX 

69, 22.97, 26.25, 32 
6 7 8 

.2 (. 2 ) .2 

.3 (.3) (.3) 


70°, 75°, 80° 


XXX 

81, 39.37 x 10^/meter 
10 12 x 10° /Ft 

.2 .2 


NOTE: Conditions in ( ) only for NACA 0008, 0012, 0015 


Boundary -layer transition strips on wings and body nose 
for all runs. 

a Range from -2° to 36°+ 

S REF = S W 

Lref = Root chord of irregular planform -C^ 

X =70 percent of Root chord -C 1 
CG J- 
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( 6 ) 


Test ARC-12-257 . Ames Research Center 12-Foot Pressure 
Tunnel. Airfoil Section Study at High Reynolds Num- 
bers and constant Mach Numbers. Limited Planform Ma- 
trix - large models (twice size of small models). 


AIRFOIL SECTION A„ = 45°; 

w 

a f = 

ON 

o 

o 

70 

o 

Ui 

O 

O 

o 

oo 

NACA 0008 



X 

X 

X 

X 

NACA 0012 



X 

X 

X 

X 

NACA 0015 



X 

X 

X 

X 

NACA 65A012 




X 

X 

X 

NACA 65]_412 
UNIT REYNOLDS 




X 

X 

X 

/" 

NUMBERS 6 . 

56, 13.12, 

16.40, 

19. 

69, 

22.97, 

26.25 X 10 /Ft 


2* 4* 

5, 


(6) 

7, 

8* 

MACH NUMBER 

.3 .3 

• 3, 

( 

.3) 

.3, 

. 3* 


NOTES: (1) NACA 65A012 and .NACA 65-^412 configurations 

run only for conditions in ( ). 

(2) Conditions noted by * run for only one con- 
figuration . 

Free transition for all runs. 

^ Range - 2° to 20° 

S REF = S W 

Lr^f = Root Chord of irregular planform -C-^ 

Xcc “ 70 percent of Root Chord 

The fact that three different test facilities and two dif- 
ferent test techniques as well as two different size models 
were used during the experimental program was considered in the 
analysis. As a consequence, there are uncertainties in the ex- 
perimental data which are reflected in the prediction methods. 

The use of different sets of reference quantities in the 
data reduction for the additional tests from those used in the 
basic planform series test required that pertinent data from 
the additional test be recalculated to be put on the same basis 
as data from the basic test. 
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TABLE 1 


GEOMETRIC CHARACTERISTICS OF BASIC WIHGS 


WHIG I 


Leading-edge sweep, deg ..... 25 

Trailing-edge sweep, deg 25 

Quarter-chord sweep, deg ..... 13.121*27 

Half-chord sweep, deg 0.00000 

Aspect ratio . 2.26500 

Taper ratio 30682 

Planfonn area, ft^ (nT) 5272U (,014898) 

Span, ft (*) 1.09279 (.33307) 

Root chord, ft (m) 73726 (.221*71) 

Tip chord, ft (m) 22768 (.06939) 

Mean aerodynamic chord, ft (n) 52733 (-16073) 

Longitudinal location of mean aerodynamic chord (x), ft (m) . . .IOU 97 (.03199) 

Spanwise location of mean aerodynamic chord (y), ft (m) 22511 (.06861) 

Airfoil sections HACA 0008, 0012 


WHIG II 


Leading-edge sweep, deg 

Trailing-edge sweep, deg 

Quarter-chord sweep, deg 

Half-chord sweep, deg 

Aspect ratio 

Taper ratio 

PI an form area, ft 2 (m2) 

Span, ft (m) 

Root chord, ft (m) 

Tip chord, ft (m) ................ 

Mean aerodynamic chord, ft(m) 

Longitudinal location of mean aerodynamic chord (x), ft (m) . . 
Spanwise location of mean aerodynamic chord (y), ft (m) . . . . 
Airfoil sections 


• • - • 35 

. . . . 20 

. . . . 23.1*6871 

. . . . 9.5li388 

2.26500 
.21*798 
.52721* (.01*898) 
1.09279 (.33307) 
.77320 (.23567) 
.19171* (.0581*1*) 
. 51*088 (. 161 * 86 ) 
.15287 (.01*659) 
.21832 (.O665I*) 

HACA 0008, Double 
Wedge 
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TABLE 1 


GEOMETRIC CHARACTERISTICS OF BASIC WIHGS 

WING III 

Leading-edge sweep, deg «... b 5 

Trailing— edge sweep, deg 15 

Quarter-chord sweep, deg . 3b. 33357 

Half-chord sweep, deg . 2Q.10b85 

Aspect ratio . 2.26500 

Taper ratio .l6bl6 

Planform area, ft 2 (s- ) . .5272b (,Ob898) 

Span, ft (m) 1.09279 (.33307) 

Root Chord, ft .(b) 82887 (.25263) 

Tip chord, ft (a) ' 13607 (.0blb7) 

Mean aerodynamic chord, ft (a) 56538 (.17232) 

Longitudinal location of mean aerodynamic chord (x), ft (m) . . .20782 (.0633b) 

Spanwise location of mean aerodynamic chord (y), ft (m) 35351 (.0633b) 

Airfoil sections 2IACA 0008 


WING IV 

Leading-edge sweep, deg 53 

Trailing-edge sweep, deg 7 

Quarter- chord sweep, deg . b 3. 9^733 

Half-chord sweep, deg 31.055b5 

Aspect ratio .. 2.26500 

Taper ratio .... 09833b 

Planform area, ft 2 (m 2 ) 5272b (.0b898) 

Span, ft (m) . 1.09279 (.33307) 

Root chord, ft (m) 87856 (.26778) 

Tip chord, ft (m) 08639 (.02633) 

Mean aerodynamic chord, ft (o) .59086 (.18009) 

Longitudinal location of mean aerodynamic chord (x), ft (m) . • .26339 (.03026) 

Spanwise location of mean aerodynamic chord (y), ft (m) 198bb (.060b8) 

Airfoil sections NACA 0008 



TABLE 1 


GEOMETRIC CHARACTERISTICS OF BASIC WIHGS 

vine V 

Leading-edge sweep, deg 60 

Tralllng-edge sweep, deg 0 

Quarter- chord sweep, deg 52.1*1091 

Half-chord sweep, deg 1*0.89615 

Aspect ratio . 2.26500 

Taper ratio .00969 

Planform area, ft 2 (m 2 ) 52721* (.01*898) 

Span, ft (m) 1.09279 (.33307) 

Root chord, ft (m) ..... . .95366 (.29128) 

Tip chord, ft (m) 00926 (.00282) 

Mean aerodynamic chord, ft (a) 63717 (. 19^21 ) 

Longitudinal location of Bean aerodynamic chord (x), ft (b) . . .31850 (.09708) 

Spanvise location of mean aerodynamic chord (y), ft (n) .... .18389 (.05605) 

Airfoil sections HACA 0008, Double 

Wedge 
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IC CHARACTERISTICS OP WING-FILLET COMBINATIONS 
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Figure 1. - Continued 














Figure 1. - Continued 





Figure 1.- Concluded 













GEOMETRIC TERMS AND EQUATIONS USED IN 
ANALYSES AND PREDICTION METHODS 


Many geometric quantities were used in the analysis and in 
development of the various prediction methods. The terminology 
and equations needed to compute specific values are presented in 
this section. 


Wing Description 


Figure 3 illustrates a typical irregular planform. To 
generate the families of irregular planforms used in this study, 
the irregular planforms were considered to be made up of a basic 
tapered planform, W, to which various fillets were added. The 
following seven geometric parameters provide sufficient informa- 
tion to allow all other planform parameters to be calculated: 

A LE’ A TE’ A F’ S W’ X W’ and (c)t/c 

\ / max 

The analyses and prediction methods consider the irregular 
planforms to be made up of an inboard panel, W,, and an outboard 
panel, W 2 . The inboard panel consists of the, fillet and that 
portion of the basic wing inboard of the intersection of the fil- 
let leading edge and the basic wing leading edge. The outboard 
panel consists of the portion of the basic wing outboard of the 
intersection of the leading edges. 


Geometric Equations 


Taper ratio of basic wing 
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Figure 3. - Wing Planform Geometry Definitions 



Nondimens ional spanwise location of break in leading edge 

sweep 

% = dy 1 /b/ 2 (2) 

Root chord of basic wing 


C R = 


/S w (tan A le + 

l 


tan A te ) 


Tip chord of basic wing 


C T X W C R 


Semispan of basic wing 

S T 


s = b/2 = 


W 


( C R + c T ) 


Span of basic wing 



Semi span of inboard panel 
dy l = 

Semispan of outboard panel 



( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 
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Chordwise distance from fillet apex to intersection of 
fillet leading edge with basic wing leading edge 

X i = C F = (dy 1 ) tan A p (9) 

Chordwise distance from apex of basic wing to leading 
edge of basic wing at the intersection with the fillet lead- 
ing edge 


C 3 = dy x tan A LE 


( 10 ) 


Chordwise distance from trailing edge of basic wing to 
trailing edge of root chord of outboard panel 

C 4 = dy 1 tan A T£ (11) 

Root chord of inboard panel 

1 - C 1 " C R + C F * C 3 < 12 > 

Root chord of outboard panel 


C 2 “ C R - < C 3 + V 


(13) 


Chordwise distance of half chord sweep of inboard panel 
across inboard panel 


C- 

0 




(14) 


Chordwise distance of half chord sweep of basic wing 
across inboard panel 


C 


6 



(15) 


Chordwise distance of quarter chord sweep of inboard 
panel across inboard panel 


C y = .75C 1 - (.75C 2 + C 4 ) 


(16) 
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Chordwise distance of quarter chord sweep of basic wing 
across inboard panel 


'8 


. 75C, 


(. 75C 0 +■ C A ) 


(17) 


Chordwise distance of maximum thickness sweep of in- 
board panel across inboard panel 


- (9 w . 


ipax 




t/c 


max 


C 2 + C 4 


(18) 


Chordwise distance of maximum thickness sweep of basic 
wing across inboard panel 


'10 


1 . 


/c 


max 


'R 


Sweep angle of quarter-chora of inboard panel 




/c. 


max 


C 2 + C 4 


(19) 


-1 


( A C/,) = tan (C 7 /dy 1 ) 

* 1 


( 20 ) 


panel 


Sweep angle of quarter-chord of basic wing or outboard 


-1 


(Ac/ 4) = tan (C g /dy 1 ) 


( 21 ) 


panel 


Sweep angle of half- chord of inboard panel 

( a C/ 2 )i “ tan" 1 (C 5 /dy 1 ) (22) 

Sweep angle of half- chord of basic wing or outboard 


-1 


(Ac/ 2) 2 = tan (C^/dy^ 


(23) 
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Sweep angle of maximum thickness line of inboard panel 

< At/c ,na* ) 1 = tan_1 (C 9 /dy l ) (24) 

Sweep of maximum thickness line of basic wing or out- 
board panel 

< At/c max> 2 = tan_1 (C i 0 /dy l ) (25) 

Average chord of inboard panel 


Cav 1 = | (C x + C 2 ) 


(26) 


Average chord of outboard panel 


: AV 2 2 (c 2 + C T ) 


(27) 


( a C/4) 


Effective quarter chord sweep of irregular planform 

( a C/ 4) ! (C AV )(dyp + ( a C/ 4) 2 (C^) (dy 2 ) (28) 


eff 


(C AV 1 ^ (dy l ) + ^ C AV 2 )(dy 2 ) 


Cosine of effective half chord sweep of irregular plan- 


form 


(Cos C/2) eff = 


Cos ( A C/2) 


l] (C AV x )(dy l 


) + 


fcos(Ac/2) 


^ C AV 2 ^ ^ dy 2^ 


< c av >< dy i> + (c Av 2 )(dy2) 


(29) 


Total planform area of irregular planform 
S TRUE = 2 ^ C AV 1 )(dy l ) + (C AV 2 )(dy 2 ) 


(30) 
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( 31 ) 


Aspect ratio of irregular planform 

2 

^RUE b / S TRUE 
Aspect ratio of basic wing 

^ " b " /S W 

Aspect ratio of inboard planform 

a = < 2d yi> 2 2d yi 

2( ' C AV^ ^ dy P C AV X 

Mean geometric chord of basic wing 



Mean geometric chord of inboard panel 


'W. 


. ^ c 
3 1 


2 1 


1 + 


W 

x+CC 2 / Cl ) 


Mean geometric chord of outboard panel 

2 


v 2 - SS 


1+ 


(c T /c 2 ) 


1+(C T /C 2 ) 

Mean geometric chord of irregular planform 


°WF = — 


(S)K) <*>l) + (%) ( C AV 2 X Jy2 ) 


(C AV 1 ) (dy l ) + < C AV 2 ><dy 2> 


(32) 


(33) 


(34) 


(35) 


(36) 


(37) 
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Spanwise location of mean geometric chord of basic wing 


y W 


1 

3 


i + 2 a w ) 
i + a w ) 


(b/2) 


(38) 


Spanwise location of mean geometric chord of inboard 

panel 



1 1 + 2<C 2 /C 1 ) 

3 r+ <c 2 /c 1 > 


(dy,) 


(39) 


Spanwise location of mean geometric chord of outboard 
panel from root of outboard panel 



X 1+2(C t /C 2 ) 

3 _r+ (c T /c 2 ) 


(dy 2 ) 


(40) 


Spanwise location of mean geometric chord of irregular 
planform 


r-t 

§ 

1 — 1 
£ 

' i-i 

) 

( dyi ) H 

'N 

(V 

CM 

X) ^ 

( 

S) 

1 1 

■) + (‘ 

' AV 2 ) 

h) 



(41) 


Chordwise location of leading edge of mean geometric 
chord of basic wing from apex of basic wing 

*W " y W 
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tanA LE 


(42) 



Chordwise distance of leading edge of mean geometric 
chord of irregular planform from apex of fillet 


*HF 


(\) (‘"'f) ( C A»J ( d >'l) + f( d Fl) f «rtp) + ( ; '» 2 K t “''LE)|( C iV 2 )( a: '2) 

Fl) f >) + f“ 2 ) W (43) 

Leading-edge radius of mean geometric chord of basic 


wing 


LER- = (100 

vJ 


(tS) 


(44) 


Leading edge radius of wing tip chord 
LER„ = (100 r/C) 

Chordwise location of moment reference point 






(45) 


(46) 


Body Description 

The models had a "minimum body” on the upper wing surface to 
house a strain-gage balance. The body consisted of a nose fair- 
ing, the balance housing and a constant crosssection aft exten- 
sion which reached to the wing trailing edge at the centerline. 
Body geometric parameters are shown in Figure 4. The nose fair- 
ing contour lines consisted of circular arcs in the longitudinal 
direction in both the vertical and horizontal planes. The cross- 
sections of the balance housing and aft extension consisted of 
flat vertical sides plus a half circle. The nose fairing cross- 
sections matched the balance housing at the point of tangency. 
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The geometric data needed for the analyses are the surface 
wetted area and the effective fineness ratio of the bodies. 


The body wetted area can be approximated by: 



Note : 


(47) 


The effective fineness ratio is defined by: 


FR = 



(48) 


Wing Wetted Area 


Inboard Panels 

S WET 1 = 4 ( C AV 1 ) ( dy l) - L N ^3 R B ) ' ( L B ‘ L n) ( 2R b) (49) 

Outboard Panels 
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EVALUATION OF EXISTING METHODS 


The initial task in this investigation was to evaluate exist- 
ing methods of predicting the aerodynamic characteristics of wings 
at low speed. Most of the existing aerodynamic prediction meth- 
ods were developed for thin sharp-edged wings of slender planform, 
for moderately thick wings of variable- sweep planform, or for mod- 
erately thick to thick wings of conventional planform. None of 
these methods was expected to apply directly to the configura- 
tions tested in the SHIPS program without additional correlation 
effort. Thus, the objective of the study of existing methods was 
to define the limits of applicability when applied to moderately 
thick to thick wings of slender irregular planform. 

While methods were sought for predicting all of the steady- 
state forces and moments in pitch, the fundamental requirement 
was to accurately predict the lift and pitching-moment behavior 
in the region of angle of attack pertinent to approach and land- 
ing of an advanced aerospace vehicle. Previous NASA studies indi- 
cated this region to be between 15 and 25 degrees of angle of at- 
tack. For that region, a significant amount of nonlinear or vor- 
tex lift occurs with many of the planforms tested. Therefore, em- 
phasis was placed on examining the ability of existing methods to 
predict nonlinear lift characteristics. 

Early theoretical analyses of the lift produced by "thin" 
wings were based on linearized potential-flow theory. These a- 
nalyses produced estimates of the slope of the lift curve (CL a ) 
evaluated near zero lift which compared well with test data over 
a limited range of angle of attack. Deviations from the linear 
extension of the lift-curve slope at higher angles of attack were 
attributed to airfoil-thickness effects and flow-separation ef- 
fects (Figure 5). In recent years, lifting-surface analyses ac- 
counting for the true surface boundary conditions have been de- 
veloped. These analyses show that the potential-f low lift (Ct ) 
is actually nonlinear with angle of attack such that the curve F 
falls below the linear estimate at all angles of attack above the 
value for zero lift (Figure 6). In addition, many wings, and par- 
ticularly those with sharp leading edges, generate leading-edge 
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flow separation in the form of bubbles (low- sweep planforms) or 
vortices (swept planforms) and tip vortices which produce lift 
above that predicted by potential- flow analyses (Figure 7). As 
the angle of attack is increased, a point is reached at which the 
bubbles or vortices burst and the lift increases at a lower rate 
and finally decreases. The angle of attack at which bursting oc- 
curs varies with wing planform and airfoil- section geometry. 

The preceding discussion indicates that a lift curve is pri- 
marily nonlinear. If a particular set of test data shows a 
linear variation over a significant range of angle of attack, it 
is apparent now that opposing nonlinear effects compensate each 
other over that range . 

The following existing prediction methods were examined: 

1. The WINSTAN nonlinear lift method for slender double- 
delta wings with sharp leading edges (AFFDL TR-66-73, 
ref. 5). 

2. The Peckham method (WINSTAN) for low-aspect-ratio 
irregular-planform wings with sharp leading edges (ref. 5) . 

3. The Peckham method as modified by Ericsson for slender 
delta-planform wings (AIAA Paper 76-19, ref. 6). 

4. The Polhamus leading-edge suction analogy for slender 
delta wings as modified by Benepe for round-leading- 
edge airfoils (GDFW ERR-FW 799, ref. 7). 

5. The empirically based Aeromodule computer procedure as 
modified by Schemensky (AFFDL- TR- 73- 144 , ref. 8). 

6. The vortex lattice lifting-surface theory computer pro- 
cedure (including effects of leading-edge suction anal- 
ogy) developed by Mendenhall et al. (NASA CR 2473, ref. 
9). 

7. The crossflow drag method for predicting nonlinear lift 
(NASA TN D -1374 , ref. 10). 

The different methods were applied either to predict the lift 
curves for a series of planforms for comparison with test data or 


49 


to attempt correlation of the test data. The results for each 
method are presented in the order listed. 


WINSTAN Empirical Method for 
Nonlinear Lift of Double-Delta Planforms 


This method is based on an empirical correlation of wind 
tunnel-test data accomplished during the WINSTAN proj ect in 1965 
(ref. 5). The method also appears in the USAF Stability and Con- 
trol DATCOM (ref. 11 ) . The calculation chart is presented in Fig- 
ure 8. Note that the method requires a value of the linear lift- 
curve slope as an input. 

As a preliminary effort, the Spencer method (ref. 12 ) was 
used to estimate Cl^ f° r a H 35 SHIPS planforms and compared with 
test data for a unit Reynolds number of 8 million per foot (eval- 
uated near zero lift). The results are presented in Figure 9. The 
Spencer method produces a satisfactory preliminary design esti- 
mate for but additional improvement is possible by applying 

a simple correction factor. The deviations of predicted values 
from test data are attributed to interference effects of the mini- 
mum bodies of the wind-tunnel models as noted by displacement of 
Cl amd near zero a, especially for the lower fineness ratio 
bodies . 

It was anticipated that the WINSTAN nonlinear method would 
predict lift coefficients higher than the test values, since it is 
known that wings with round-leading-edge airfoils produce less 
nonlinear lift than do wings with sharp-leading-edge airfoils. 

THE WINSTAN nonlinear lift prediction method was applied to 
two series of SHIPS planforms. One series consisted of the 80- 
degree leading-edge-sweep fillet and outboard-panel sweeps of 
25 to 60 degrees. The other series consisted of an outboard- 
panel sweep of 25 degrees with various fillet sweeps from 80 to 
35 degrees. 

Figure 10 presents predicted nonlinear and linear lift 
variations with angle of attack for the 80-degree-sweep-fillet 
series of planforms. Also shown is the range of experimental 
data for this series of planforms at 16 degrees of angle of at- 
tack. The experimental data falls between the linear and non- 
linear estimates for all the outboard-panel sweeps. 
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All of these planforms have true aspect ratios which are less 
than 1.65. It is apparent that the basic flow field is dominated 
by the leading-edge vortices generated by the highly swept fillets; 
thus, one must expect only a small variation of lift with outboard- 
panel sweep at a constant angle of attack. Figure 11 illustrates 
this fact. Note that the round- leading- edge airfoils produce 
about 50 to 55% of the vortex lift that is estimated by the WIN- 
STAN method for sharp leading edges. 

The second series in which the fillet sweep is varied pro- 
duces considerably larger variations in the lift curves. Figure 
12 presents estimates for this series of planforms. The results 
of the WINSTAN nonlinear method are shown for all the planforms. 
Note that the planforms with fillet sweeps of 35, 45, and 55 de- 
grees are outside the bounds of the data base for the WINSTAN non- 
linear correlation, but the estimates appear quite reasonable. 

Also shown are linear estimates of lift and the nonlinear 
potential-flow lift for the three of the planforms. The ranges of SHIPS 
test data are indicated for 12 and 16 degrees of angle of attack. 
For fillet sweeps from 35 to 70 degrees, the test data fall below 
the linear estimates but above the estimates. Figure 13 

shows the variation of lift coefficient with fillet sweep for a 
constant angle of attack of 16 degrees to illustrate this re- 
sult more clearly. It is apparent that little vortex lift is 
produced for fillet sweeps of 60 degrees or less with the 25- 
degree-swefep outboard panel, which suggests a boundary condition 
exists with respect to the fillet contribution to nonlinear lift. 

It is of interest to plot the ratio of vortex lift produced 
by the SHIPS wings to the vortex lift predicted by the WINSTAN 
nonlinear method for sharp- leading-edge wings. Figure 14 shows 
that when a strong vortex is produced by a high-sweep fillet, 
there is little variation in the vortex lift at 16 degrees of 
angle of attack with sweep of the outboard panel. The variation 
of the vortex-lift ratio with fillet sweep clearly shows that the 
vortex lift is small up to 60 degrees of fillet sweep and then 
increases rapidly as fillet sweep is increased to 75 degrees. 

There is apparently a plateau slightly above 75 degrees. 

As expected, the WINSTAN lift prediction method for sharp- 
leading-edge wings significantly overpredicts the lift generated 
by the SHIPS models. However there are two possible ways to use 
the method as a basis for a new method. One way is to derive 
families of plots similar to those of Figure 14 for various 
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angles of attack, fillet sweep angles, outboard- panel sweep an- 
gles, and Reynolds numbers and then to seek analytical expres- 
sions to curve-fit the data. A second way would be to use the 
WINSTAN nonlinear-method correlation parameters to define a cor- 
relation of the SHIPS data for irregular planforms with round- 
leading-edge airfoils. Reynolds number effects would require a 
separate correlation chart. 


The Peckham Method for Low- Aspect-Ratio 
Irregular-Planform Wings with Sharp Leading Edges 


The Peckham method as applied to irregular-planf orm wings 
was also investigated during the WINSTAN program (ref. 5). This 
method, which is most appropriate for slender wings, consists of 
a plot of the correlation parameter Cl/(s//)^ against angle of 
attack, as shown in Figure 15. The parameter s is the wing 
semispan and l is the overall length parallel to the plane of 
s ymme t ry . 

The Peckham method was applied to correlate SHIPS test data 
for three series of planforms, as shown in Figures 16, 17, and 18. 
The spread of data in these curves is systematic with either 
inboard- or outboard-panel sweep. It is apparent that the term 
(s / l)^ is too powerful for the round-leading-edge wings. 

As an alternate, a weaker relationship, (s//) 4 , was at- 
tempted. Figure 19 presents the correlation achieved with the 
parameter Cl/(s//)^. Note that below (X = 15°, data for all the 
fillet sweeps from 35 to 80 degrees can be represented by a single 
curve, whereas above Ot = 15° there is still a systematic spread 
in the data. This spread is obviously caused by the variation in 
vortex breakdown effects with sweep angle; and, thus, it should be 
possible to gain some insight about correlating the effects from 
the body of literature covering vortex breakdown. The modified 
Peckham method was thus considered a candidate for additional 
evaluation . 
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The Peckham Method as Modified by Ericsson 
for Slender Delta-Planform Wings 


Ericsson’s version of the Peckham correlation method for 
slender delta wings (ref. 6) uses the parameters 


C L /(A/4) 3/2 = f (a/e LE) 

where A is the aspect ratio, ot the angle of attack, and ^LE is 
the complement of the leading-edge sweep angle. Note that for a 
delta wing, A/ 4 = s/l = tan ^LE • By making the assumption that 
\an& LE ~ for slender wings, Ericsson divides each side of the 

Peckham correlation by^LE- 

Several attempts were made to define an effective for ir- 

regular planforms that would correlate the test data using Eric- 
sson’s approach. None were successful. The best of the several 
that were tried used area-weighted cosine ^ LE to define the #lE 
effective. The results are presented in Figure 20 for six of the 
SHIPS planforms. It is apparent that dividing the angle of at- 
tack by the effective #lE cutises the data to spread rather than 
collapse because the outboard-panel- sweep contribution varies 
rapidly as the sweep changes. This approach was dropped from 
further consideration. 


The Polhamus Leading-Edge Suction Analogy 
as Modified by Benepe for Round-Leading-Edge Airfoils 


This method, which was originally developed for round- leading- 
edge delta wings, is summarized in Figure 21. The method makes 
use of the leading-edge suction parameter ”R, " which is defined 
by the equation: 


„„„ - C L tan “ + CD min ' C » 
R ~ — r o 

C T tana - L 
7TA 
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The numerator in this definition represents the equivalent amount 
of suction actually present while the denominator represents the 
amount of suction theoretically possible for a wing of the same 
aspect ratio but having an elliptic loading. In the past, "R" 
has been successfully used in analyzing and predicting drag-due- 
to-lift in the lift coefficient range up to polar break. 


In the present investigation, it had been proposed to utilize 
the suction ratio concept to define the nonlinear lift of round 
leading-edge wings by combining the "R" concept with Polhamus 
leading-edge suction analogy concept (ref. 12) which is applicable 
to sharp leading-edge wings. This has successfully been done for 
low to moderate aspect ratio delta wings (ref. 7) . 


The lift prediction approach is represented by the following 
equations : 


C 


L 



+ 5 


C L 

VOR 


SLE 



where C T is the nonlinear potential flow lift which is de- 
bp 

fined by: 


C T = CL sin O' cos a 
bp b# 


and 8 


C L 


IS 


VOR 


SLE 


the vortex lift produced by an equivalent plan- 


form having an infinitely thin sharp leading-edge airfoil. Ac- 
cording to Polhamus leading-edge suction analogy concept, the 
leading-edge suction calculated by potential flow theory actually 
produces a normal force on a sharp thin highly swept wing. 


In the present evaluation of existing methods, the 5. 


VOR 

term is obtained by subtracting C from values of CL predicte 

Lp L 


from a correlation of lift of sharp-edged double delta irregular 
planforms developed prior to the revelation of Polhamus leading- 
edge suction analogy concept of vortex lift during the WINSTAN 
project studies (ref. 5). Thus, 


6 
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For convenience, the C, term in the nonlinear potential lift is 

Let 

evaluated by Spencer's method (ref. 12) for irregular p lan forms : 

277A 1 

C T 


ot 


2 + 


prr 

A 

2 

J 

(cosA c/2)eff 



57.3 per degree 


For wings having round leading-edge airfoils, it is assumed 
that the difference between full leading-edge suction and the ac- 
tual leading-edge suction is converted to vortex lift. Thus, the 
term ji-"R (O0"J is applied to to obtain the vortex 


VOR 


SLE 


lift increments. 


This method was applied to analyze data for several SHIPS 
planforms. The variations of "R" with angle of attack calculated 
from the test data are presented in Figures 22-24. These varia- 
tions are typical of the results that can be expected from the 
SHIPS data. Figure 22 presents data for Basic Wing V at four dif- 
ferent unit Reynolds numbers; it illustrates the fact that a sig- 
nificant Reynolds number effect occurs throughout the angle-of- 
attack range. Figure 23 presents data for two basic planforms 
and two irregular planforms -- all obtained at the same unit 
Reynolds number, 3.98 x 10 6 per foot. Note the significant dif- 
ference in the shape of the variation of "R" with a. for Basrc 
Wings I and V. The variation shown for Wing I is representative 
of upper- surface flow separation developing on a low-sweep wing 
and the loss of lift associated with the flow separation. The 
variation shown for Wing V, on the other hand, is representative 
of the generation of a vortex-type flow separation, which gener- 
ates additional lift above that for potential flow. The varia- 
tion for the combinations of Wing I with a 60-degree fillet shows 
the general characteristics of that for Basic Wing I; and, thus, 
one should expect little vortex lift to be generated at high 
angle of attack. The variation shown for Wing I with the 80- 
degree fillet^ also exhibits the dominant effect of the low-sweep 
outboard panel; but, since a lower level of "R" exists at low to 
moderate angles of attack, some vortex lift should be present 
on that wing. 

Figure 24 presents the variations of "R" with Of for the Basic 
Wing V and Wing V with an 80-degree fillet. The more rapid de- 
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crease in "R" at low to moderate angle of attack is indicative of 
a stronger vortex flow than that occurring on the Basic Wing V, 
and one would expect a large nonlinear lift contribution for this 
irregular planform. 

The test variations of "R" with a. were then used to predict 
the lift curves for four of the planforms. The linear lift- curve 
slope was predicted by Spencer's method. The nonlinear 
potential-f low curves were generated by use of the linear lift-curve 
slope and the appropriate geometric equation for correcting the 
linear lift to nonlinear lift. Sharp-edged vortex-lift incre- 
ments were obtained by estimating the total nonlinear lift for an 
equivalent sharp-edged thin wing by use of the WINSTAN method and 
subtracting the nonlinear potent ial- flow values. The effects of 
the round leading edges were thus computed by applying l-"R(ce)" 
values to the sharp-edged vortex-lift increments. The results 
are presented in Figures 25-28. 

Figure 25 is for Basic Wing V. Predictions are compared 
with test results for two values of unit Reynolds number. Both 
predictions are slightly higher than test values; however, the 

correlation of Spencer's C T prediction with test data shown in 

a 

Figure 9 indicate, d the predicted value to be about 6 % high. A 67 0 
reduction in the estimates for 20 degrees of angle of attack in- 
dicates good correlation between the test data and the estimates. 
Figures 26 and 27 show the predictions for Wing I with 60- and 
80-degree fillets. The predictions are obviously poor above the 
angle of attack for stall of the low-sweep outboard wing panel. 

It thus would be necessary to modify the analysis /correlation 
method to account for the fact that loss of leading-edge suction 
on the low-sweep outboard panels does not necessarily contribute 
additional lift. 

Figure 28 presents the result for Wing V with an 80-degree 
fillet. In this case, the prediction compares well with test. 

While th£re were some difficulties associated with this 
method, it was attractive from the standpoint that defining the 
variations of "R" with a with the planform geometry parameters, 
and with leading-edge-radius Reynolds number was also considered 
the best approach for predicting the drag-due-to-lift . 
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At this point, several of the SHIPS models were altered to 
obtain force and moment data with the outborad wing panels removed 
so that the leading-edge suction parameter variations with ot could 
be determined for use in estimating the fillet contributions. The 
results of Test 8TPT-780 (inboard panel alone) were initially 
referenced to the total planform area of the irregular planforms 
from which they were made. The present analysis required that the 
lift and drag values be re-referenced to the actual planform areas . 

Figures 29 and 30 show typical results for two different 
fillet sweep angles. Figure 20 presents variations of the leading- 
edge suction Ratio - "R" for the inboard panel alone, the basic 
wing and the combined irregular planform for the 80- degree fillet 
and Wing I with the NACA 0008 airfoil section. On the basis of 
previous analyses for delta wings, the variation for the inboard 
panel was expected to have a generally similar shape to that of 
the irregular planform, but lower values of "R" at each angle of 
attack. The actual variation was entirely different from what 
was expected, starting with a low value at low angle of attack 
and increasing with angle of attack. Figure 30 presents a similar 
comparison for the 65-degree fillet and Wing I, and the results 
show the same type of effects. After considerable contemplation 
of these results, it was concluded that the tip vortex formed 
along the side edge of the inboard panel models interacts favor- 
ably with the leading-edge vortex to enhance the leading-edge 
suction as angle of attack is increased. 

This concept was verified by analyzing test data for delta 
and cropped delta wings having thin, biconvex airfoils with sharp 
leading edges (ref. 14). The results are presented in Figure 31. 
Note that for the cropped wing, the values of leading-edge suction 
ratio are initially higher than those for the delta wing and 
increase at higher angles of attack while values for the delta 
wing decrease throughout the angle-of -attack range. 

Further examples of the leading-edge suction ratio data 
and lift prediction comparisons are presented in Figures 32 
and 33, respectively, for the irregular planform consisting of 
Wing III with 70-degree sweep fillet. Results are shown for 
both the NACA 0008 airfoil section and the thin modified double 
wedge airfoil section models. 
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Figure 32 shows that the variation of "R" with Of is signifi- 
cantly different for the wings having an NACA 0008 airfoil and a 
sharp leading-edge modified double wedge airfoil. Note that the 
forward facing slopes of the sharp leading-edge airfoil produce 
significant suction values at low angle of attack. The value of 
"R" decreases rapidly with increasing a for the sharp wings. At 
a = 14.3 degrees a kink was also noted from the basic data (not 
presented) . The naviation for the round-leading-edge airfoil is 
much different. "R" remains at relatively high values up to 11.4 
degrees a. . Above that a , a distinct change occurs in the slope 
of the curve. These changes, noted for either airfoil, are 
important clues to a change in flow field which causes the amount 
of suction to be less than the value 1 - "R(a)"J. 

Figure 33 presents the buildup of lift predictions and com- 
parisons with test data for the two wings previously discussed. 

The curves labeled WINSTAN S.L.E., SPENCER LINEAR and C T repre- 

Lp 

sent elements of the prediction method that are dependent only 
on the planform and -Mach number. The curves labeled SLE and RLE 
have the term |l- n R(ot) , j applied to determine the vortex lift as 
explained earlier. 

For the sharp-edged airfoil, the agreement between prediction 
and test is good up to 14.3 degrees when a change in the flow 
field apparently occurs (probably an outboard panel stall). For 
the round- leading-edge airfoil, the agreement between prediction 
and test is good up to 11.4 degrees where again a change in flow 
field occurs although the effect on the experimental lift curve 
is not as apparent as it is for the sharp-leading-edge airfoil. 

Working the method in reverse order, the suction ratio val- 
ues required to make the predicted lift match the test data were 
calculated for these two cases and the results are shown in Fig- 
ure 34. In essence, the difference between the dashed curves 
and the test data curves of "R" represents the amount of "lost 
suction" that is not converted to vortex lift. 

In general, the nonlinear method based on the leading-edge 
suction analogy overpredicted the lift at high angles of attack 
for the round- leading-edge airfoil irregular planform wings. 

There is a trend toward better agreement over a larger range of 
angles of attack as the fillet sweep and outboard panel sweep 
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increase. The best agreement occurs, as expected, for Wing V 
with the 80-degree sweep fillet (Fig. 40). 

Reynolds number effects were evaluated for the irregular 
planform consisting of Wing I with 80 degree fillet. Both the 
NACA 0008 and NACA 0012 airfoil data were analyzed. For the 
NACA 0008 airfoil, the agreement between prediction and test 
data improved with decreasing Reynolds number because the test 
vortex lift increased faster than the predicted vortex lift. 

For the NACA 0012 airfoil, the agreement between prediction and 
test improved as the Reynolds number increased because the pre- 
dicted vortex lift decreased faster than the test vortex lift. 

(The specific data substantiating these statements are not pre- 
sented, however the reader is referred to the appendix where the 
incremental effects of Reynolds number on the finalized lift 
correlation parameter are presented. The effects are highly 
configuration dependent.) 

It is apparent that for this irregular planform, the effects 
of Reynolds number are threefold. First, the angle of attack for 
"stall" of the outboard wing panel is increased as the leading- 
edge-radius Reynolds number is increased. Second, the strength 
of the leading-edge vortex flow produced on the inboard panel 
is reduced substantially as the leading-edge radius Reynolds 
number is increased. Third, the interaction between the inboard 
and outboard wing panel flow changes drastically between the two 
extremes of leading-edge-radius Reynolds number tested. 

The obvious difference in inboard panel/outboard panel flow 
interactions is that at low Reynolds number, the outboard panel 
having the NACA 0008 airfoil produces a leading-edge flow 
separation, while at high Reynolds number, the outboard panel 
having the NACA 0012 airfoil produces a trailing-edge flow separa- 
tion. It is likely that the wing having NACA 0008 airfoils pro- 
duces leading-edge flow separation on the outboard panel at all 
unit Reynolds numbers tested, whereas it is possible that mixed 
leading-edge and trailing-edge flow separation occurs for the 
NACA 0012 airfoil at the lowest unit Reynolds number tested. 

From this evaluation of applying the leading-edge suction 
analogy to predict nonlinear lift of the SHIPS planforms, it 
was concluded that a significant amount of resources and time 
would be required to make the approach a useful preliminary de- 
sign tool. 


59 


Aeromodule Computer Procedure 


The Aeromodule computer procedure was developed to provide a 
rapid means of assessing the lift, drag and pitching moment char- 
acteristics of aircraft configurations, from low speed to super- 
sonic speeds. The basis for the calculations is an extensive 
set of empirical and semi-empirical methods developed from a data 
base encompassing fixed wing and variable sweep wing configura- 
tions. The computer program has been revised and updated several 
times since the initial development in 1969 as the data base ex- 
panded and imprpved prediction methods were developed. A version 
of the program (General Dynamics Procedure Code R1T) developed 
for the U.S. Air Force (ref. 8) was used in the present investi- 
gation. 

Aeromodule predictions were obtained for most of the SHIPS 
planforms during this evaluation. Comparisons between Aeromodule 
predictions of lift , drag, and pitching moments with the SHIPS 
wind-tunnel data obtained at 26.25 million per meter unit Rey- 
nolds number are presented for a mini-matrix of 9 configurations 
shown in Table 3. These results adequately illustrate some ma- 
jor points to be made from the Aeromodule prediction evaluation 
study. 


TABLE 3 MINI-MATRIX OF CONFIGURATIONS PRESENTED 


Fillet Sweep \ 

Wing Leading-Edge Sweep 

- Deg. ! 

Deg T 

25 

45 

60 | 

j 

None 

Figure 35 

Figure 36 

! 

Figure 37 j 

80 ; 

Figure 38 

Figure 39 

Figure 40 

65 


Figure 43 

55 

1 

Figure 42 


1 

35 J 

Figure 41 

— 

| 

i 

i 
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The overall evaluation for the Aeromodule prediction methods 
indicated that, while there were sometimes significant differences 
between predicted and test values, the methodologies contained 
the basic elements that were needed to provide more accurate pre- 
dictions and could be adapted through data correlations to a- 
chieve the necessary accuracy. In general, the discrepancies 
could be explained by the fact that the SHIPS wings and test con- 
ditions fall outside the data base that was available when the 
empirical methods were developed. 

Comparisons for basic Wings I, III, and V are presented in 
Figures 35, 36, and 37, respectively. These results were infor- 
mative for they showed that the empirical factors inherent in 
the computer procedure were not completely satisfactory for the 
basic planforms. This is not surprising, since the SHIPS basic 
wings are generally lower in aspect ratio or have thicker air- 
foils than wings which supplied the data base. In addition, the 
26.25 million per meter unit Reynolds number was higher than pre- 
viously available in any significant parametric study. 

The differences between the predicted lift curves and drag 
curves and the test data were not distressing from a methodology 
development standpoint. The SHIPS data base could provide ade- 
quate information to define new empirical factors. The pitching- 
moment predictions were actually encouraging even though the dis- 
crepancies are quite large. 

The fact that the general nonlinear character of the pitching- 
moment curves also occurred in the predictions was a signifi- 
cant point. Obviously, some alteration to the methodology was 
required, but again, it could only be a matter of changing the 
empirical factors to better account for stall progression on the 
wings with angle of attack. This assumption is borne out by 
the results shown (Fig. 38, 38, and 40). 

Figures 38, 39, and 40 present the comparisons between pre- 
diction and test for Wings I, III and V with 80 -degree fillets. 
For all of these irregular SHIPS planforms, the Aeromodule com- 
puter program logic selected the WINSTAN nonlinear lift predic- 
tion method for sharp leading-edge wings, which does not in- 
clude the effects of round leading edges or Reynolds number. The 
predicted lift curves tend toward better agreement with test 
data as the outer wing panel sweep increases to 60 degrees be- 
cause the round- leading-edge and Reynolds number effects decrease 
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There is an improvement in the agreement between the pitch- 
ing moments and test data corresponding to the improvement in 
agreement of the lift predictions. 

The drag-polar predictions for the wings with 80-degree 
sweep fillets all produced higher drag values than test at mod- 
erate lift coefficients. These results were also due to the lack 
of accounting for round- leading-edge and Reynolds number effects. 

The third set of 3 configurations for which comparisons are 
presented between Aeromodule predictions and SHIPS test data con- 
sists of Wing Iwith the 35-degree sweep fillet, Wing III with 55- 
degree sweep fillet, and Wing V with the 65-degree sweep fillet. 
The comparisons are shown in Figures 41, 42 and 43, respectively. 
For the first two cases, the computer program logic initially 
selected the WINSTAN nonlinear lift method for double-delta plan- 
forms, but then attempted to match the USAF DATCOM low-aspect- 
ratio CL max method results and determined that the CL max value 
predicted by that method was lower than the lift predicted by po- 
tential flow. The program flow then shifted to the WINSTAN non- 
linear lift method for cranked wings and the DATCOM high-aspect- 
ratio C;L max met h°d despite the fact that the configuration is a 
low-aspect-ratio wing. In this case, the procedure computed a 
CL max but did not attempt to fair from the computed lift curve 
to The user must do that by hand. 

The corresponding moment curves in Figures 41 and 42 reflect 

the use of the Ct values above the angle of attack at which 

J max 

the computed lift values exceed . The pitching-rnoment pre- 

diction methodology assumes a center of pressure location for a fully 
separated flow at lift values above CL max - A better prediction 

of the angle of attack at which Ct occurs could improve the 
° L max r 

calculated pitching moments. 

It was not at all obvious from the comparisons shown in 
Figure 43 (foi; Wing V with 65-degree sweep fillet) what computa- 
tional path was used by the computer procedure, but it appeared 
to be the same one as occurred for the configurations presented 
in Figures 41 and 42. Unfortunately, the specified C^ range for 
the series of predictions only extended to a value of 1.0, so the 
pitching-moment curve does not give a clue. 


62 



The comparisons between predicted and test drag polars pre- 
sented in Figures 41, 42 and 43 again show that the predicted 
drag is too high in the moderate C L range and too low in the 
higher C L range. Improved accounting of Reynolds number effects 
on the lift coefficient for initial flow separation and the lift 
coefficient for drag break could lead to significant improvement 
in both these areas. 

The longitudinal stability derivative dC m /dC L is compared 
in Figure 44 with data evaluated in the low lift coefficient 
range for Wing I combined with several different fillets. The 
methodology is that of Paniszczyn from the WINSTAN study (ref. 5) 
which also appears in the USAF Stability and Control DATCOM (ref. 

11) in the section on wing-body pit chin g -moment predictions for 
double- delta , cranked and curved planforms. The discrepancies 
which occur are attributable in part to discrepancies in pre- 
dicting the lift curve slope. The only drawback to the method 
is the fact that it requires generation of geometry for a ficti- 
tious outboard panel. 

Comparisons between predicted and test values of the minimum 
drag coefficient, > are presented in Figure 45. In general 

the discrepancies between predicted values and test values were 
much larger than was anticipated from previous experience. The 
variations of the experimental data with fillet sweep show some 
significant deviations from smooth curves. Some of this type of 
discrepancy was found to be attributable to the effects of the dif- 
ferent forebody shapes on the measured base pressures. In addition, 
the inherent lack of sensitivity at low angle of attack of a wind- 
tunnel balance when designed to measure large forces at high angle 
of attack contributed to the "scatter. " Investigation as to why 
the predicted values were so much higher than expected revealed 
that the geometry calculations contained in the Aeromodule program 
could not orouerly account for the fact that the body existed only 
on the upper surface, so the calculated wetted areas for the body 
contributions to minimum drag were too large. When the predicted 
values were corrected to account for the actual wetted areas, the 
agreement between prediction and test values was much better and 
within the expected accuracy. 

In summary, an extensive evaluation of the Aeromodule pre- 
diction methods was completed, and the conclusion was that the 
methods contain most of the necessary elements to produce accu- 
rate predictions for the SHIPS planforms. The discrepancies 
that did occur are primarily due to an inadequate data base 
which the SHIPS test data could remedy. The Aeromodule methods 
were used as guidelines during the data correlation efforts. 
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Lifting- Surface Theory of Mendenhall, et al. , 
Including Effect of Leading-Edge Suction Analogy 


This method (ref. 9) was included in a preliminary study be- 
cause it offers a convenient way to obtain theoretical values of 
the nonlinear lift contribution for sharp leading edges and also 
because the amount of leading-edge suction converted to lift can 
be arbitrarily applied for each wing panel of a two-panel irregular- 
planform wing. Analysis of Wing I with a 60-degree fillet 
was tried, first, by using full- leading-edge suction analogy and, 
then, by arbitrarily assigning various amounts of leading-edge 
suction to the inboard and outboard panels at each angle of at- 
tack. Guidance in selecting the applied amounts was obtained 
from the plot of R(a') for the configuration. The results are 
presented in Figure 46. 

The plot shows the variation of , the WINSTAN sharp- 
leading-edge prediction, the Mendenhall predictions for full and ar- 
bitrary amounts of suction assigned, and the test data. The Men- 
denhall result for full suction conversion to vortex lift agrees 
with the WINSTAN prediction. The result with partial suction 
conversion agrees better than when the test value of suction was 
applied to the WINSTAN increment. The reason is that at ct = 15° 
and above the outboard-panel suction ratio was set to zero since 
lift data for Basic Wing I indicated wing stall at about a = 15°. 

The next step was to apply the same technique as had been 
done in the WINSTAN study for cranked wings of moderate thickness, 
that is, to assume that at 15 degrees of angle of attack the 
potential-f low contribution had ceased to grow and the only addi- 
tional lift was that caused by the vortex flow on the inboard 
panel. This result is shown by the symbol x in Figure 46. The 
agreement with test data is excellent. The choice of 15 degrees 
was arbitrary, in this case; what is needed is a valid correlation 
of a for as a function of planform parameters, airfoil 

suction parameters, and Reynolds number. The basis for such cor- 
relations is contained in the Aeromodule methods, but data perti- 
nent to the SHIPS planforms would be helpful to refine the 
approach. 
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The Crossflow Drag Method 



The crossflow drag method for predicting nonlinear lift 
(ref. 10) assumes that the total lift is obtained by a relation 
of the form 

C L ’ C L tf (C ' ) + C D C (-573-) 

where the term C n is the planform drag coefficient evaluated at 

C 

a = 90 at high Reynolds number and a units are degrees. 

Although this method is useful for very preliminary estimates 
of the nonlinear lift, the C~ term is relatively insensitive 

C 

to wide differences in planform shape for low- aspect-ratio wings. 
Values in the literature vary from about 1.15 to 1.30 for wings 
with round-leading-edge airfoils. In addition, the available 
data are not systematic with planform shape, so it would be diffi- 
cult to apply the method with sufficient accuracy to differen- 
tiate between the effects for the various SHIPS planforms . 
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Linear Prediction 


‘ Test Data 



Figure 5. Early Analysis Approach 



Figure 6. Recent Potential Flow (Lifting- Surface Theory) Analysis 


Test Data 



Figure 7. Analysis Including Effect of Vortex Flow 




Figure 10. Effect of Outboard-Panel Sweep on Lift 
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Figure 11. Effect of Outboard-Panel Sweep on Lift at a * 16° 
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Figure 12. WINSTAN Nonlinear Predictions Wing I Planforms 
(Sharp Thin Airfoils) 
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Figure 13. Effect of Fillet Sweep on Lift of 
Wing Planforms at a * 16 
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Figure 14. Effect of Round Leading Edges on Vortex Lift 
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Figure 15. Correlation of Lift Curves of Gothic and 
Ogee Planforms 
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Figure 16. Effect of Outboard Panel on Peckham-Type 
C L Correlation (A f * 80°) 
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Figure 23. Variations of Suction Ratio with Angle of Attack Figure 24. Variations of Sucticn Ratio With Angle of Attack 

for Various Planforms at Unit Reynolds Number of Showing Effect of Fillet 
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Figjre 27. Prediction of Nonlinear Lift Using Test Values Figure 28. Prediction of Nonlinear Lift UEing Test Values 

of Suction Ratio for Wing I with 80° Fillet of Suction Ratio for Wing V with 80° Fillet 
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Wing III With 70° Fillet 
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Fieure 34. Modified Variations of Leading-Edge- Suction Ratio 
Wing III With 70° Fillet 



Figure 35, Comparisons of Aeromodule Predictions With Test 
Data - Configuration 25 25.0008 































Figure 39. Comparisons of Aeromodule Predictions With Test 
Data - Configuration 80 45.0008 








Figure 40. Comparisons of Aeromodule Predictions With Test 
Data - Configuration 80 60.0008 













Figure 42. Comparisons of Aeromodule Predictions With Test 
Data - Configuration 55 45.0008 






Figure 43. Comparisons of Aeromodule Predictions With Test 
Data - Configuration 65 60.0008 
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Figure 44. Comparison of Aeromodule Longitudinal Stability 

Derivative With Test Data for Wing I With Various 
Fillets 
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Figure 46. Results of Analysis Using Mendenhall Lifting 
Surface Plus Suction Analogy Computer Proce- 
dure and WINSTAN Methods 



DATA CORRELATIONS 


The goals of the data correlation efforts were twofold: 

(1) to provide a means of removing the deficiencies 
and limitations of existing prediction methods 
or developing new empirical methods, and (2) to 
generate simplified design guides for evaluating 
new configuration concepts. 

The specific objective of the task was to obtain combina- 
tions of planform geometric parameters, airfoil section param- 
eters and flow parameters that would correlate the SHIPS lift, 
drag, and moment data. 

The evaluations of existing prediction methods indicated 
that new or revised correlations were needed for the following 
elements of the aerodynamic characteristics: 

(1) Second order (nonlinear) lift 

(2) Drag due to lift 

(3) Aerodynamic-center location 

(A) High angle -of -attack pitching moment characteristics 
(pitch-up boundaries) . 

The approach used was to obtain the basic effects of plan- 
form parameters from the test data obtained at 26.25 million per 
meter unit Reynolds number with the NACA 0008 airfoils and then 
examine the effects of Reynolds number and airfoil section para- 
meters. The efforts will be described in the order listed above. 

Nonlinear Lift 

The analysis of existing methods indicated that accurate pre- 
diction of nonlinear lift was perhaps the most essential ingre- 
dient in the development of an improved prediction methodology. 
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Therefore, a considerable amount of effort was devoted to investi 
gating different correlation approaches for lift. 

A plot showing the spread of the lift coefficient data 
for all 35 SHIPS planforms with NACA 0008 airfoils is pre- 
sented in Figure 47. The coefficients are referenced to the to- 
tal planform area. What is desired is to find correlation para- 
meters that will coalesce these data into either a single line or 
a family of lines. From previous work it was known that plotting 
data for constant values of fillet sweep would reduce scatter con 
siderably. Figure 48 which presents the data for SHIPS planforms 
with 80 degree fillets illustrates this fact. Note that the re- 
maining scatter is systematic with outboard panel sweep. A fur- 
ther reduction in scatter for angles of attack below 20 degrees 
was obtained by dividing the test lift coefficient (C^) by the 
lift curve slope per degree angle of attack C^q- (estimated bv the 
Spencer method) as illustrated in Figure 49. Figure 50 presents 
the same correlation parameter applied to the lift data for the 
Basic Wings. The reduction in scatter is excellent below 16 de- 
grees angle of attack. 

The WINSTAN correlation parameter for lift of double delta 
planforms was examined next. The parameter combines C^/Cl^ with 
two irregular planform geometric parameters: (1) A^ the aspect 

ratio of the inboard panel and (2 ) t ?b, the spanwise location of 
the break in leading-edge sweep as a fraction of semispan. The 
parameter is 



In this case, the lift-curve slope is per radian. 


Figures 51 through 58 present correlations of SHIPS lift 
data for each of the fillet sweeps starting with the 80- degree 
sweep fillet wing combinations. In each case, the correlation 
produces a satisfactory collapse of the data up to 16 degrees 
angle of attack and a reasonably systematic spread of the data 
above that angle, which is dependent on outboard panel sweep. 


The data were then further correlated using the second cor- 
relation parameter from the WINSTAN double-delta lift-prediction 
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method. Curves of/ C, 



as functions of tanXp were estab- 


C T V 

LtOC 

lished at constant 'angles 6f attack as shown in Figure 59. The 
curves for angles of attack from 0 to 16 degrees apply to all 
outer panel sweeps. The curves above 16 degrees apply only 
for = 45 degrees. 


Further effort to account for other outboard panel sweeps 
was held in abeyance until other correlations of lift were a- 
chieved. It was obvious, however, that the WINSTAN double-delta 
method was a prime candidate for developing a SHIPS lift predic- 
tion technique. 


The next lift correlation technique investigated was the 
modified Peckham method in which the parameter Cl/(s//)% i- s 
plotted against angle of attack. Figure 60 is a plot 
showing the spread of data for all 35 SHIPS planforms. The col- 
lapse of data below 16 degrees angle of attack is remarkable con- 
sidering the simplicity of this correlation. The spread of data 
above 16 degrees reflects, of course, the differences in stall 
progression among the many wing planforms. Figures 61, 62 and 63 
present individual correlations for the basic wings and for the 
fillet wing combinations with 65-degree and 80-degree fillets. 

The data collapse up to 16 degrees angle of attack is not 
significantly better than that of the overall correlation. 

The data above 16 degrees do indicate definite families of 
curves for the various outboard panel sweeps. 

A further modification of the Peckham method was also evalu- 
ated. The previously described correlation parameter was divided 
by the lift-curve slope, per degree, as predicted by the Spencer 
method. The new parameter 



is plotted against angle of attack in Figures 64 through 72 for 
the basic wings and then for wing fillet combinations for constant 
values of fillet sweep. The data collapse for angles of attack 
below 16 degrees is much improved compared to using only 
except for fillet sweeps of 65 and 70 degrees where the (s/i)^ 
improvement is less significant. In addition, correlation curves 
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I 

4 


'f 


representing these results have different initial slopes. A plot 
similar to the WINSTAN double-delta method presented earlier would 
be needed to have a viable prediction technique. Although the cor- 
relations do not appear to be quite as good as those of the WIN- 
STAN method, this approach did offer an alternative to the WIN- 
STAN method. 

The many different approaches presented thus far were aimed 
specifically at trying to develop correlations that could evolve 
into a prediction method. The next two analyses were intended 
primarily to evolve design guidelines. 

First, the SHIPS lift data were plotted using the para- 
meter 



SPENCER 


in which C T is the linear estimate of lift variation with 

angle of attack. Figures 73 through 81 present the plotted cal- 
culations for the basic wings first and then the fillet-wing com- 
binations for constant fillet sweeps. The primary use of these 
curves would be to establish the upper limit of angle of attack 
for quasi-linear lift for each configuration. Rather than pro- 
vide arbitrary boundaries, it was decided to allow the designer 
to exercise some judgment as to how much nonlinearity would be 
acceptable for a given situation. The design guide would con- 
sist of faired curves grouped like the test data plots. 


The final correlation parameter examined is the ratio of 
test lift coefficient to the nonlinear potential flow lift co- 
efficient , C T /C T 

L Lip 

where 

2 

C T = C, cos a sin Of and C_ is the Spencer value 
L p L a La 
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per degree. Typical variations of this parameter with angle of 
attack are shown in Figures 82 and 83. The results for 80 de- 
gree sweep fillets (Figure 83) show large contributions due to 
vortex lift, whereas the results for the basic wings exhibit lit' 
tie vortex lift except on the 60 -degree sweep wing. 


Drag Due to Lift 

The basic approach to developing correlations for drag due 
to lift was the analysis of the leading-edge- suction ratio, M R" , 
as a function of angle of attack, planform parameters, airfoil 
section parameters, and Reynolds number. The basis of the ap- 
proach was the fact that for low to moderate angles of attack, 
correlation of leading-edge-suction ratio with planform para- 
meters and the effective leading- edge radius Reynolds number had 
been achieved and verified for both swept and irregular planforms 
during the WINSTAN investigation (Ref. 5). 

As an initial step in the analysis , . values of the suction ra- 
tio were calculated for all runs of test ARC 12-086-1 and plotted 
versus angle of attack. Figure 84 presents a typical family of 
suction ratio data for the SHIPS planforms with 75 degree fillet 
sweeps obtained at unit Reynolds number of 26.25 million per 
meter (8 million per foot) and illustrates a fact noted in most 
of the plotted data. There are six distinct regions to the varia- 
tions with angle of attack. At low angles of attack the varia- 
tions are erratic because the analysis is extremely sensitive to 
the value selected for Cd^^ and the test values of lift and drag 
coefficients are not sufficiently accurate to define specific 
data trends. At angles above about 4 degrees, a plateau 
value is apparent which has been used in previous work to define 
the basic drag polar shape for attached flow. For analysis pur- 
poses, the lower limit of the plateau region was defined as OL ^ 
and the upper limit as ol ^ • Next there is a slight decrease in 
the slope of the variations which is apparently caused by the 
initial development of flow separation on the wing. While this 
region is probably a curve, for analysis purposes it is con- 
sidered to be linear, the upper limit of this region is defined 
as The flow separation then develops more rapidly causing 

a marked decrease in the value of the suction ratio with angle of 
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attack, which again appears to be nearly linear. The upper limit 
of this region is defined as OC^. Next, there is a reduction in 
the rate of decrease of the suction ratio with angle of attack as 
the spread of flow separation slows. The upper limit of this re- 
gion, which also apparently marks the beginning of a region in 
which the upper surface flow is fully separated, is defined as ct^. 

As an initial step in the analysis the values of the angles 
of attack bounding each of the regions were determined from the 
test data obtained on all the planforms at unit Reynolds number 
of 26 .25 /million per meter. Typical results are presented in 
Figure 85 which shows the variations of and 

with fillet sweep angle for the Wing I planforms. 

Further examination of the variations of suction ratio with 
angle of attack for the wings with NACA 0008 airfoils suggested 
that for prediction purposes the slopes of the variations in a 
given region could be considered to be the same at each of the 
Reynolds numbers for a particular planform. Thus it might be 
possible to establish a prediction method which could account 
for Reynolds number effects on the angle -of -attack boundaries. 

It was first necessary to establish a correlation of the 
plateau values with planform geometry and Reynolds number. Fur- 
ther consideration of the plotted data led to the following 
approach: 

(1) A correlation of suction ratio would be obtained 
for the basic wing planforms. 

(2) The effect of fillet sweep would be accounted for 
by a correction term. 

Two correlations of the basic wing data were achieved. 

The first correlation uses the effective leading-edge radius 
Reynolds number in which the velocity and radius are both taken 
normal to the, swept leading edge. The results are compared in 
Figure 86 with the band of data presented by W. P. Henderson in 
NASA TN D-3584 (Ref. 15) for wings of 4 to 6 percent thickness 
ratio and boundary- layer transition fixed by the Bras low tech- 
nique (Ref. 16). The SHIPS data agree quite well with the earlier 
data, but there- is slightly more scatter. The reason for the 
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scatter is most likely the fact that the balance axial- force sen- 
sitivity is lower than those used for the data presented in 
TN D-3584, because the balance had to be capable of handling loads 
at high angles of attack, high values of dynamic pressure and a- 
cross the speed range. 

The second correlation used the WINSTAN correlation para- 
meter (Ref 5) . 

E ^ER C0TA LE }J l ~ H C ° S A LE 

where: R^^ is t ^ ie leading-edge radius Reynolds number based on 

the streamwise velocity and radius measured at the mean geometric 
chord. The other two terms provide an empirical fit for sweep 
and Mach number effects. The original WINSTAN correlation also 
used another parameter aA to further account for planform 

effects. cosA t „ 


The results of applying this approach to the SHIPS basic 
wing test data are presented in Figure 87 along with the original 
WINSTAN correlation curves for AX values of 0 and 1 which 


cosA t f 

bound the SHIPS planform values. Tne SHIPS data which were ob- 
tained with fixed transition form a different correlation band 
than the original WINSTAN data which were obtained with free 
transition. In addition, the planform correlation does not ap- 
pear to be appropriate for the SHIPS data. 


Figure 88 presents a comparison of correlated values of suc- 
tion ratio with test data assuming the long dash short dash line 
in Figure 87 represents a correlation curve for the effects of 
Reynolds number and leading-edge sweep, "R" BW > for the basic wings. 

Using the premise that the data correlation curve, "R'' y , 
shown in Figure 87 should form the basis for evaluation of tne 
effects of Reynolds number on the leading-edge suction ratio at 
low lift coefficients for the irregular planforms , it was possible 
to establish families of straight line "curves" (shown in Figure 
89) with fillet sweep angle and unit Reynolds number as the para- 
meters . 
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The families of curves presented in Figure 89 can be repre- 
sented by another family of curves shown in Figure 90 which present 
the curve slopes A"R" as functions of outer panel leading- edge 

sweep angle, Ay, for parametric values of leading- edge radius 


Reynolds number based on streamwise velocity and radius measured 
at the mean geometric chord of the basic wing. Thus, the lead- 
ing-edge suction ratio at low lift coefficients can be predicted 
using a relatively simple equation: 


Uplt = tin" 

1 0<oc<a 2 


'R' 


BW 


~ A”R" 

(Ap “Ay) 


(Ap ”Ay) 


where "R" w is obtained from Figure 87 and the term in brackets 
is obtained from Figure 90. For simplicity the "plateau value" 
of "R" is assumed to apply from zero angle of attack to the boun- 
dary a 2 - 


With a satisfactory methodology in hand for predicting the 
suction ratio in the plateau region, attention was then turned 
to refining the evaluation of the angle-of -attack boundaries and 
slopes for the other regions. Despite several iterations between 
selections of the boundaries and slopes.no satisfactory correla- 
tion could be achieved when using the boundaries directly, except 
forc^* The curves defining Of ^ are presented in Figure 91. How- 
ever, it was finally noted that incremental values of angle of 
attack between each boundary formed reasonable sets of curves 
for regions 2 and 3 as shown in Figure 92 and 93. The incre- 
mental values for region 4 were somewhat erratic but definable 
as shown in Figure 94. 


It therefore appeared that a workable prediction method 
for the variation of suction ratio with angle-of-attack was 
achievable if a suitable correlation could be found for the ef- 
fect of Reynolds number on the angle-of-attack boundaries. That 
aspect of the .methodology was completed in the following way. 


First, plots were made of the values of the ot boundaries 
as functions of fillet sweep with parametric variations of unit 
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Reynolds number. Figure 95 is a typical example. After much 
thought, it was decided to attempt to correlate the boundaries 
for the basic wings as functions of Reynolds number and then ap- 
ply a ratio term to the increment in a for the basic wings to ac- 
count for the effect of fillet planforms . 

A clue to a possible method of correlating the data for the 
basic wings was obtained from a paper by Chappell of the Royal 
Aero Society's Engineering Services Data Unit (Ref. 17). He 
had correlated the angle of attack for initial separation, a s , 
of swept, tapered wings by plotting ttgCosApE against the leading- 
edge radius Reynolds number evaluated normal to the wing at the 
tip. Figure 96 shows f° r each of the basic wings correlated 
using Chappell’s approach and comparisons of his results for air- 
foil thickness ratios of .08 and .10. The results were encourag- 
ing.. 

Then a form of the Q function was tried in which the stream- 
wise leading-edge radius Reynolds number at the tip was substi- 
tuted for the value based on the mean geometric chord. The re- 
sults are presented in Figure 97. This approach was tried with 
all the ot boundaries. There was enough consistency to the re- 
sulting plots to attempt to pursue this approach. 

Since a correction term was sought rather than the absolute 
values, the individual values for each basic wing sweep were 
shifted along the ordinate axis to form a reasonably smooth en- 
velope curve encompassing the data obtained at the maximum unit 
Reynolds number of 26.25 million per meter (8.0 million per 
foot) . Figures 98 through 101 show the results for each angle- 
of-attack boundary. The correction terms A(S#( ^cosAy) due to 
Reynolds number effects are obtained from Figures 98 through 101 
by first calculating values of£p at a reference Reynolds number 
and the Reynolds number appropriate to the conditions to be 
evaluated, then subtracting the value of (8a( )COsA y) ref from 
the value (5a( ^cosA^) as illustrated in Figure 98. 

The reference value of Reynolds number is evaluated assuming 
a unit Reynolds number of 26.25 million per meter (8.0 million 
per foot) and a tip leading- edge radius corresponding to a wing 
span of 0.33307 meters (1.09279 feet). If the reference value of 
Pp is less than the value of Pp for the condition to be evaluated, 
then A(8a( ^cosA y) will have a positive value. 
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The next step was to obtain the effect of the irregular 
planforms by forming ratios of the increments in angle of attack 
due to Reynolds number at the various values of fillet sweep to 
the increments measured for each basic wing planform. The re- 
sults are shown in Figure 102. These ratios are applied to the 
increments produced by the previous step to obtain the final re- 
sult for the effect of Reynolds number on the angle-of-attack 
boundaries for each region. 
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Aerodynamic -Center Location 


The initial step in the pitching -moment analysis was to sur- 
vey the - C.^ curves for each configuration for the data ob- 


tained from test 12-086 at unit Reynolds number of 26.247 million 
per meter (8.0 million per foot). In this initial survey, 3 dis- 
tinct quasilinear regions were observed for many of the wings: 

A small region near zero lift, a ’’primary -slope” region generally 
falling between lift coefficients of 0.1 and 0.5, and a region at 
high lift coefficients. 


The aerodynamic -center locations (referenced to the mean 
geometric chord of each planform) for each region are presented 
in Figures 103 through 107. The angle of attack and lift coeffi- 
cient envelopes for the primary slope region are shown in Figures 
108 and 109. In general, the variations of a.c. location with 
fillet sweep in the low-lift and "primary-slope" regions are 
smooth curves, whereas the variations for the high-lift region are 
irregular for fillet sweeps in the neighborhood of 60 to 70 de- 
grees. The envelope curves for the upper and lower limits of 
the "primary-slope" region also shows distinct irregularities in 
!:he region of fillet sweep from 60 to 70 degrees. It is apparent 
that when a significant amount of vortex flow is present on the 
fillet there is a definite reduction in the shift of a.c. loca- 
tion between the primary-slope region and the high angle -of -at- 
tack region. In a few cases, the shift in a.c. location is 
actually destabilizing at higher lift coefficients. 


Emphasis was then placed on attempting to correlate the a.c. 
location in the "primary-slope" region. One approach was to 
modify the Panisczcyn method of predicting the a.c. location (de- 
veloped in the WINSTAN investigation, Ref. 5) by applying a mul- 
tiplying factor k to the lift prediction for the inboard panel. 

A parametric evaluation was conducted for a range of k values 
from 0 to 2. The results are presented in Figures 110 through 
114. The comparisons with test data indicated that weighting 
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factors (k) in the range from .9 to 1.1 would give good correla- 
tion. 


In seeking an alternative to the Panisczcyn method which 
might be more simple to apply, the a.c. locations in the "pri- 
mary-slope" region were analyzed and plotted in several different 
formats as presented in Figures 115 through 121. None of these 
formats provided a significant collapse of the data, but in each 
case the data form families of curves for parametric values of 
either fillet sweep or outboard panel sweep. Eventually, Figures 
120 and 121 were curve fitted to produce simple and useful pre- 
diction techniques. 

The effects of Reynolds number on the a.c. location in the 
low and high-lift-coefficient regions and the "primary- slope" 
regions are shown in Figures 122 through 126. The effects were 
small in both the low-lift and "primary- slope" regions for all 
configurations, but are quite pronounced for some of the confi- 
gurations at high lift coefficients. In general, an increase in 
unit Reynolds number produces an aft shift in a.c. location. 
Figure 125 for example shows that for Wing IV with 65 degree 
fillet an increase in unit Reynolds number from 19.67 million 
per meter (6.0 million per foot) to 26.25 million per million(8.0 
million per foot) changes a pitch-up trend to a pitch-down trend. 

As a consequence, efforts to analyze pitch -up/pitch-down 
characteristics were initially concentrated on the data obtained 
at 26.25 million per meter unit Reynolds number. 


High Angle-of-Attack Pit chin g -Moment Characteristics 


As an aid to describing the pitching -moment characterisitcs , 
the - C L curves were reduced in size and assembled on one page 

for all 35 planforms (Figure 127). At low angles of attack, the 
change in the static stability derivative dC^/dC^ is in the nega- 
tive direction with increasing values of either fillet sweep or 
outboard panel sweep. At high angles of attack, the changes in 
dC m /d^ are in the opposite direction. All of the planforms with 
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fillet sweep of 60 degrees or less and outboard panel sweep of 
53 degrees or less exhibit a definite pitch-down tendency at high 
angles of attack. As the fillet sweep or the outboard panel 
sweep is increased, the trend at high angle of attack is toward 
less stability until eventually a definite pitch-up tendency oc- 
curs. The moment reference point is 0.25cyp in Figure 127. 

The pitch-up/pitch-down tendencies are summarized in Figure 
128 which is a matrix chart having the same arrangement as Figure 
12 7 in terms of the placement of the data for each combination 
of fillet sweep and outboard panel sweep. Superimposed on the 
chart are arbitrary boundaries which separate the regions of de- 
finite pitch-down and definite pitch-up from the other configura- 
tions. In the region between these boundaries, the data does not 
define whether pitch-up or pitch-down exists. 


the a BREAK and C L 


Some attempts were made to quantify and correlate the shift 

in a.c. location from the "primary slope" region to higher 

lift coefficients and the angle of attack and lift coefficient 

for the break point (upper limit of the primary slope region) of 

the pitching moment curves. A typical example is shown in Figure 

129 in which Aa.c. BDrAiri Of " d id t a and C T are plotted against 

BREAK BREAK L B REAK 

the variable ^tan (Ap -Ay ) . This approach produced groupings of 

data which showed some promise, so effort 

BREAK 

was concentrated on the CL analysis and several correlation 

BREAK 

parameters were tried and the most promising results (shown in 
Figure 130) used the parameters ^^TRUE and tan " Ay )• 

BREAK 

Figure 131 presents the same correlation approach after the 
CL points had been reviewed and considerable smoothing was 

BREAK 

applied to the raw results. The smoothed data show two distinct- 
ly different trends. Planforms having either Wing I or Wing II 
outer panels ,show a decreasing value of C, /A-,^— as 

L BREAK TRUE 

tan (Ap - Ay ) increases, whereas planforms having Wing IV and 


Wing V outer panels show an increasing value of 


/A 


BREAK 


TRUE * 
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Planforms having Wing III outer panels apparently are on the 
borderline between the two trends . 

Although this correlation was informative, it did not appear 
to be easily translated into a prediction method. Shortly 
thereafter it was noted that the angle-of -attack boundaries ob- 
tained in the suction-ratio analysis corresponded closely the 
break points occurring the plots of pitching -moment variation 
with angle of attack for many of the planforms. Figures 132 
and 133 show typical examples of fitting the experimental pitch- 
ing moment curves with linear segments between the angle -of - 
attack boundaries established from the suction-ratio analysis. 

This approach was adopted and it was found that two addi- 
tional angle-of -attack boundaries were needed for the pitching- 
moment analysis: which corresponded to the lower limit of 

the plateau region of the suction -ratio analysis and OL ^ which 
occurred at high angle of attack for a few configurations. Thus, 
it appeared possible to produce prediction methods for lift, drag, 
and pitching-moment variations with angle of attack which would 
produce consistent results. 

The analysis of pitching-moment variations with angle of 
attack also brought to mind an alternate approach to developing 
a design guide for the high angle-of-attack pitching behavior. 

The author remembered a comment in a British paper on the use 
of irregular planforms to minimize the aerodynamic- center travel 
which occurs between subsonic and supersonic flight. The comment 
warned that along with the favorable reduction in a.c. travel 

came a tendency toward pitch-up at high angles of attack at low 
speed. 

This author interpreted that comment to mean that a corre- 
lation between a.c. location at low lift coefficients and high 
angle-of-attack pitchup might exist. Such a correlation was at- 
tempted, and the results presented in Figure 134 show that it is 
possible to define a boundary between pitch-up and pitch-down be- 
havior as a function of outboard panel sweep and aerodynamic- 
center location in terms of percent mean geometric chord of the 
irregular planforms. 
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One of the goals of the overall SHIPS program is to provide 
a guide for using irregular planforms yielding desired improved 
Cl, reduced transonic a.c. shift, while avoiding pitch-up. Figures 
128 and 134 provide one element of the overall design guide appro- 
priate to the landing condition. Note that the design guides do 
not rule out any configuration, they merely point out that cer- 
tain configurations must be evaluated in more detail than others. 
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Figure 47. Scatter Plot Showing Spread of Lift Data for all 
35 SHIPS Planforms 
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Figure 48. Plot Showing Spread of Lift Data for Constant Fillet 
Sweep 









ure 49 . Examp 





























Figure 51. Correlation of SHIPS Test Data Using WINSTAN Corre- 
lation Parameter for Nonlinear Lift of Double-Delta 
Wings - Wings with 80 Fillet Sweep 




igure 52. Correlation of SHIPS Test Data Using WINSTAN Corre- 
lation Parameter for Nonlinear Lift of Double-Delta 
Wings - Wings With 75°- Sweep Fillets 
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Wings - Wings With 70 - Sweep Fillets 



Wings - Wings With 65 - Sweep Fillets 





Figure 55. Correlation of SHIPS Test Data Using WINSTAN Corre- 
lation Parameter for Nonlinear Lift of Double-Delta 
Wings - Wings With 60 - Sweep Fillets 
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Figure 56. Correlation of SHIPS Test Data Using WINSTAN Corre- 
lation Parameter for Nonlinear Lift of Double-Delta 
Wings - Wings with 55 - Sweep Fillets and Basic Wing 
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Figure 57. Correlation of SHIPS Test Data Using WINSTAN Corre- 
lation Parameter for Nonlinear Lift of Double-Delta 
Wings - Wings With 45 - Sweep Fillets 
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Figure 58. Correlation of SHIPS Test Data Using WINSTAN Corre- 
lation Parameter for Nonlinear Lift of Double-Delta 
Wings - Wings With 35° Sweep Fillets and Basic Wing I 
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Figure 59. SHIPS Lift Data Correlation Using WINSTAN Method 
for Double-Delta Wings (Preliminary) 




Figure 60. Scatter Plot for Modified Peckham Lift Correlation 
Parameter for All 35 SHIPS Planforms 
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Figure 61. Modified Peckham Correlation of Lift Data - Basic 
Wings 
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igure 62. Modified Peckham Correlation of Lift Data 
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Figure 64. Data Collapse With Further Modification of Peckham 
Lift Correlation Parameter - Basic Wings 
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Figure 68. Data Collapse With Further Modification of Peckham 
Lift Correlation Parameter - A„ = 65 
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Ratio of Actual Lift to Predicted Linear Lift of SHIPS 
Planforms With NACA 0008 Airfoils - Basic Wings 









[PS Planforms with NACA 0008 Airfoils 
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Figure 76. Ratio of Test Lift to Predicted Lineai Lift of 
SHIPS Planforms with NACA 0008 Airfoils - A„ = 
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Figure 77. Ratio of Test Lift to Predicted Linear Lift of 
SHIPS Planforms with NACA 0008 Airfoils - A„ = 
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Figure 82. Ratio of Test Lift to Nonlinear Potential Flow Esti- 
mate for Basic SHIPS Planforms 
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CONFIGURATION 
75 25 .0008 



Figure 84. Variations of Leading-Edge-Suction Ratio with Angle 
of Attack for Irregular Planforms with 75° Fillet 
Sweep 
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Figure 86. Correlation of Suction Ratio "R" Using Effective 
Leading-Edge Radius Reynolds Number 
















Figure 87. Correlation of Suction Ratio Using WINSTAN Q Parameter 












Figure 88. Comparison of Correlation Curve Values of Suction 
Ratio with SHIPS Basic Wing Test Data 









Figure 90. Calculation Chart for Incremental Effect of Fillet 
Sweep on Suction Ratio in Plateau Region Between 
a = 0° and <X 
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Figure 96. Correlation of Reynolds Number Effects on a 2 for 

Basic Wings Using Chappell's Correlation Parameter 
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Figure 98. Envelope Correlation Curve for Change in °f 
Basic Wings Due to Change in Reynolds Number 
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Figure 99. Envelope Correlation Curve for Change in a 3 1 
Basic Wings Due to Change in Reynolds Number 
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SHIPS Wind-Tunnel Data 


M = .3, R N = 8.0xl0°/Ft 



Figure 108. Angle-of-Attack Envelope for Primary Aerodynamic- 
Center Location 



SHIPS Wind-Tunnel Data 
M = .3, R N = 8.0xl0 6 /Ft 



Figure 109. Lift-Coefficient Envelope for Primary Aerodynamic- 
Center Location 
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Figure 113. Aerodynamic Center Prediction - Paniszczyn Method 
Modified by K Factor Applied to Lift of Inboard 
Panel - Wing IV Planforms 
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Figure 1X4 


Aerodynamic- Center Prediction - Paniszczyn Method 
Modified by K Factor Applied to Lift of Inobard 
Planforms - Wing V Planforms 
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SHIPS Wind-Tunnel Data 
M = .3, Rn = 8 . 0xl0 6 /Ft 
Primary Slope 
(Cl~ ■ 1 to .5) 



Figure 115. Aerodynamic -Center Location Referenced to Mean 
Geometric Chord of Each Irregular Planform 
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SHIPS Wind-Tunnel Data 

M = .3 
35° Wing 



Figure 123. Effect of Reynolds Number on Aerodynamic-Center 

Location for Various Lift Regions - Wing II Plan- 
forms 
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SHIPS Wind-Tunnel Data 


M = .3 
60° Wing 



Figure 126. Effect of Reynolds Number on Aero- 
dynamic-Center Location for Various 
Lift Regions - Wing V Planforms 
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Figure 127. Curve Shape Schematics 
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M : '■■■' NACA 0008 Airfoils 


Probably_will Pitch 
Y_ Up at Higher CX — 

Unstable Break at 


-(deg) 


May Pitch Up or Down 
- — i — at Higher Ct 


Stable Break 


Figure 128. Pitch Up/Pitch Down Tendencies as 
Functions of Wing and Fillet Sweep 
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Figure 131. Correlation of Lift Coefficients for Upper Limit 
of Primary Slope Region Based on Smoothed Wind- 
Tunnel Data 
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Figure 133. Example of Variation of Pitching Moment with Angle 
of Attack Approximated by Linear Segments Between 
Angle-of'Attack Boundaries Defined by Suction-Ratio 
Analysis - Irregular Planform with High Fillet Sweep 











Figure 134. Pitch-Up/Pitch-Down Design Guide Based on Low 
Angle of Attack Aerodynamic -Center Location 
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PREDICTION METHOD DEVELOPMENT 


The previously described efforts produced several approaches 
for developing a set of prediction methods for lift, drag, and 
pitching-moment characteristics, but additional work remained. 

As discussed earlier, it was intended to devise a basic set of 
methods from the data obtained at unit Reynolds number of 26.25 
million per meter (8.0 million per foot). This section describes 
how selected approaches for lift, drag, and pitching moment (in 
that order) were developed into the "basic" prediction technique. 

One guideline in the method development process was the fact 
that eventually prediction methods would be required for tran- 
sonic and supersonic speeds and compatibility among the methods 
would be desirable. A second guideline was the desire to comput- 
erize the final methods. A third guideline quite obviously was 
the need to obtain sufficient accuracy from the predictions to be 
able to represent the differences in the aerodynamic characteris- 
tics that occurred between changes in the irregular planforms. 

The methods selected for development or modification were : 

(1) The WINSTAN nonlinear lift correlation 

(2) The Aeromodule prediction method for minimum drag 

(3) The leading-edge suction ratio correlation approach 
for drag due to lift 

(4) The quasilinear variation of pitching moment with- 
in regions of angle of attack defined from the 
drag -due -to -lift prediction 

The development efforts related to predicting the lift 
characteristics are presented first, followed in order by those 
for the drag characteristics and the pitching -moment characteris- 
tics . 
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Lift Characteristics Prediction 


While many approaches for correlating the lift characteris- 
tics of the SHIPS planforms were examined and several showed pro- 
mise, the WINSTAN lift correlation approach shown previously in 
Figure 59 was selected for further development because it had 
produced the best collapse of data at angles of attack up to 16 
degrees and contained parameters which would allow extension to 
higher Mach numbers. In addition, it was apparent the spread of 
data above 16 degrees angle of attack was sufficiently syste- 
matic in terms of the values of outboard panel leading-edge 
sweep for each value of fillet sweep that the chances of finding 
a means of accounting for the effects were good. 


The approach taken was to redefine the basic correlation 
for angles of attack above 16 degrees on the basis of data for 
the outboard panel sweep equal to 25 degrees. These incremen- 
tal values of the lift correlation parameter a / Cl Ai \ were 

determined for each value of outboard panel Icl# 7/b / 
sweep. Plots were made of the incremental values as functions 
of the parameter /StanAp for various fixed angles of attack 
from 16 degrees to 26 degrees. The results were somewhat ir- 
regular for fillet sweeps below 65 degrees and angles of at- 
tack above 20 degrees, but it proved possible to smooth the 
baseline curves in a way (Figure 135) which also produced 
smoother variations of the incremental values with /3tanAp. 
It was also noted that, if the smoothed incremental values at 
22 degrees angle of attack (Figure 136) were used as a basis, 
the variations of the incremental values with angle of attack 
at fixed values of /QtanAp could be reproduced to good accuracy 
by a correlation curve, f^, (Figure 137) which represents the 
stall progression of the family of 35 wing planforms. 

The final results provided the following prediction equa- 
tion : 


C 


L 


a 


'b 


% 




+ f 


l A 
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where the terms in the brackets are read from correlation charts 
as follows : 



A 




f 


1 


from Figure 135 


from Figure 136 


from Figure 13 7 


is the aspect ratio of the inboard panel 


V)-. is the non-dimensional spanwise location of the 

' D 


break in the leading-edge sweep. 


The term f L ^1 \ is a function of the parameters /JtanA-r, 

v* V 


and the angle of attack,# . The term ii 





is a function 


of 0tanAp andA^ that is applied only for angles of attack 

greater than 16 degrees and accounts for the effects of outboard 
panel sweep. The factor f^ varies with angle of attack above 16 

degrees to reproduce the incremental effect on stall progression 
when applied with the incremental effects caused by variations 
with outboard panel sweep. 

The lift curve slope, C, , is the incompressible value cal- 

L a 

culated using Spencer s semi- empirical expression: 
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2 7T A 


per radian 



where : 

A is the "true" aspect ratio of the irregular planform, 

and 


cos Ac/2 


eff 


cosA c/ 2 Wi * S 1 + cosA c/2 W2 * S 2 
S 1 + S 2 


The lift values calculated by this method inherently in- 
clude interference effects of the "minimum fuselages" or balance 
housings to which the various wings were attached. While the 
effects are not strictly definable, the basic result is to intro- 
duce a slight negative camber to each wing which is manifested 
in a slightly negative value of lift at zero angle of attack. 

The method also has the limitations of being applicable strictly 
to wings with NACA 0008 airfoil sections and Reynolds numbers 
appropriate to the maximum values tested during ARC Test 12-086-1. 
However, it will be shown later that the method will provide a 
reasonable estimate of full-scale lift characteristics for the 
NACA 0008 airfoils. 


Drag Characteristics Prediction 


The approach used to develop a prediction method for the 
drag characteristics of the irregular planforms consisted of 
adapting an existing methodology for minimum drag of complete 
aircraft configurations (Reference 8) and finalizing a drag- due 
to- lift methodology based on the leading-edge- suction ratio 
correlation efforts discussed earlier in this report. In both 
cases it is possible to account for Reynold’s number effects. 
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Minimum Drag 

The minimum drag of the irregular planform wings is com- 
posed of drag items that are "assumed" to be independent of lift 
such as friction, form, interference, base, roughness and protu- 
berance drag contributions. The methods used to determine each 
of the minimum drag contributions are discussed in the following 
subsections . 


Friction, Form, and Interference Drag 


A large part of the subsonic minimum drag is comprised of 
the sum of friction, form, and interference drag of all the con- 
figuration components. The drag of each component is computed 
as 



where is the compressible flat-plate skin-friction coefficient, 
S wet is the component wetted area, and FF and IF are the compo- 
nent form and interference factors. 

For emphasis, the specific configuration components con- 
sidered are illustrated in the sketch below. The pertinent geo- 
metric equations are presented in an earlier section 




Hidden area 
not used in 
wetted area 
calculation 
(Body side 
and wing 
upper surface 
only) 
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Friction Drag 

The flat-plate, compressible, turbulent, skin-friction coef- 
ficient is determined from the general equation given in Refer- 
ence 18 


- ¥[ Cf i (X • * 2 ) 


(3-2) 


where and are functions of the free-stream Mach number and 

wa ll temperature. The incompressible skin- friction coefficient, 
r , is evaluated at the equivalent Reynolds number, R^. F ? . 
f ± L 

White and Christoph (Reference 18) developed expressions for the 
transformation functions F^ and along with a more accurate 

explicit equation, based on Prandtl/Schlichting type relations, 
for computing the incompressible, turbulent, flat-plate friction 
coefficient (C f ) with the following results: 
i 

F, = f 1 f 1 

f 2 = t 14 *f 


0.430 


r-\ _ T> \*- • 36 

log 10 

For an adiabatic wall condition, t and f are given by 


1 ^ log 10 


t ='Too/T 


aw 


•[ 


l + r^ Moo 2 


-1 


f = 1 4- 0.044 r Mqo t 
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Using a recovery factor r = 0.89 and a viscosity power-law expo- 
nent n = 0.67, recommended in Reference 18, results in the follow- 
ing expression for C^: 


where 



(3-3) 


f = 1 + 0.03916 M m 2 . t 


The Reynolds number, R*. , is based either on component length or 

L 

an admissible surface roughness, whichever produces a smaller 
value of Reynolds number, as follows: 



minimum 


(Rjg/ ft) 


K. 


(L/K) 


L 

1.0489 


(3-4) 


where 


R^/ ft is determined from standard atmospheric tables 
or is input. 

L is the characteristic length of the component. 

K is the admissible surface roughness and is an 
input quantity. 


and, 

K x = 37.587 + 4.615M + 2.949M 2 + 4.132M 3 
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For mixed laminar- turbulent flow, transition location is 
specified for the upper and lower surfaces of the wing. For 
the laminar portion of the flow, the Blasius skin- friction rela 
tion 


C f = 1.328/ 



(3-5) 


Laminar 


is used up to the transition 


where C f /C f = (1 + 0. 1256M 2 )' - 12 
i 

point. At the transition point, X^, the laminar momentum thick- 
ness is matched by an interative process to a turbulent momentum 
thickness, which begins some fictitious distance, AX, ahead of 
transition. The skin- friction coefficient for the turbulent part 
of the flow is calculated from Equation 3-3, where the Reynolds 
number is calculated from 


^ = (AX + L - X r ) 


(V ft) 


(3-6) 


The value of with transition is finally given by 


C 


f 




Turb . 


(3-7) 


Calculated values of versus R„ are presented in Figures 138 

L 

through 143 for mixed laminar-turbulent flow. 


Form Factors 

The component form factors, FF, account for the increased 
skin friction caused by the supervelocities of the flow over 
the body or surface and the boundary- layer separation at the 
trailing edge. The form factor for the "body" component is 
computed as 

FF * 1 + 60/FR 3 + 0.0025 • FR (3-8) 
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where 


Component Length 
^Width x Height 

For "nacelle" components, the form factor is given by 

FF = 1 + 0. 35/FR (3-9) 

Equations 3-8 and 3-9 were obtained from the Convair Aerospace 
Handbook (Reference 19) and also appear in the DATCOM (Refer- 
ence 11) . 

The airfoil form factors depend upon airfoil type and stream- 
wise thickness ratio. For 6-series airfoils, the form factor is 
given by 

FF = 1 + 1. 44 (t/c) + 2 (t/c) 2 (3-10) 

For 4-digit airfoils, the form factor is given by 

FF = 1 + 1. 68(t/c) + 3(t/c) 2 (3-11) 

For biconvex airfoils, the form factor is given by 

FF = 1 + 1. 2 (t/c) + 100(t/c) 4 (3-12) 

And for supercritical airfoils, the form factor is given by 
FF = 1 + K.C. , + 1. 44(t/c) + 2 (t/c) 2 (3-13) 

The factor K-. C, , in Equation 3-13 is an empirical relationship 
which shifts the 6-series form-factor equation to account for 
the increased supervelocities caused by the supercritical-section 
design camber C.^. The factor (derived from experimental 

data) is shown plotted in Figure 144 as a function of the Mach 
number relative to the wing Mach critical. Equations 3-10 and 
3-11 were obtained from informal discussions with NASA/LRC per- 
sonnel; Equation 3-12 appears in both the DATCOM and the Convair 
Handbook. 
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Interference Factors 

The component interference factors, IF, account for the mu- 
tual interference between components. For the fuselage, the in- 
terference factor is given by 

IF = %-B (3 ’ U > 


where ^ is shown plotted in Figure 145 as a function of fuse- 
lage Reynolds number and Mach. For other bodies such as stores, 
canopies, landing-gear fairings, and engine nacelles, the inter- 
ference factor would be an input factor based on experimental 
experiences with similar configurations. The Convair Aerospace 
Handbook (Reference 19) recommends using 

IF = 1.0 for nacelles and stores mounted out of 
the local velocity field of the wing . 

IF = 1.25 for stores mounted symmetrically on 
the wing tip. 

IF = 1.3 for nacelles and stores if mounted in 
moderate proximity of the wing. 

IF = 1.5 for nacelles and stores mounted flush 
to the wing or fuselage. 

The interference factor for the main wing is computed as 

IF = ’ ^W-B (3-15) 

where R^_g is the wing-body interference factor presented in 

Equation 3-14, and R^g is the lifting surface interference 

factor presented in Figure 146. For supercritical wings the 
wing interference factor is set equal to one. Other airfoil 
surfaces such as horizontal or vertical tails use an interfer- 
ence factor determined by 

IF = R^g • H f (3-16) 
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where is the hinge factor obtained from input (use = 1.0 

for an all movable surface, 1.1 if the surface has a flap for 
control) . The factors R^_g and R^g are plotted in Reference 5 
and also appear in the DATCOM. 


Rase Drag 


Data presented in Reference 20 were used to establish equa- 
tions from which the base drag of bodies could be determined. 

The trends of these data show three different phases: (1) a 

gradual rise of at transonic speeds up to M = 1, (2) a 

Base 

relatively constant drag level supersonically up to about M = 
1.8, and (3) a steadily decreasing value of drag above M = 1.8. 
The resulting empirical equations are given as 


i 


8 S 

(0.1 + 0.1222M ) Base , M<1 

3 Ref 


'D 


Base 


0. 2222S„ / S „ 1.0<M<1.8 

Base Ref 


^Base^Ref ^ 3 * 15 +M >■ M>1 * 8 


(3-19) 


Miscellaneous Drag Items 

In the preliminary design stage of aircraft drag estimation, 
the drag due to surface irregularities such as gaps and mis- 
matches, fasteners, small protuberances, and leakage due to 
pressurization are estimated by adding a miscellaneous drag in- 
crement which is some percentage of the total friction, form, 
and interference drags. The miscellaneous drag varies between 
5 and 20 percent of the total friction, form, and interference 
drags for typical aircraft. 
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Drag Due To Lift. 


The approach selected for development of a drag-due- to- lift 
prediction method is based on the concept of the leading-edge 
suction ratio. The various correlations developed have already 
been discussed in some detail. The task that remained was to 
formalize a calculation procedure using forms of the correla- 
tions suitable for calculation and define the specific equa- 
tions to be used. As a consequence, some of the charts pre- 
sented with the following summary of the prediction method are 
repeats of previous figures. 

The basic equation for drag due-to-lift, , is: 

C n = C T tana - "R(a) M 

J_i 


is the drag due-to-lift with zero leading edge 
suction 

is the drag due-to-lift with full leading edge 
suction 

is the ratio of leading-edge suction derived 
from the test data to the full leading-edge 
suction value as a function of angle of attack. 

is the "true" aspect ratio of the irregular 
planform . 

The variation of M R(OD " with CX typically consists of 5 
distinct quasi-linear regions as illustrated in Figure 147. 

Region 1 is the plateau region for which M R(Q0" is a con- 
stant for a specific aircraft geometry and test condition, but 
is a function of leading-edge radius Reynolds number, Mach 


where : 

tana 



7T A 

M R(a) M 

A 


C T tana - L 

7TA 
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number, and leading-edge sweep of both the basic wing planform 
and the fillet planform. 


The value of "R^" is obtained from the equation: 


IV 


f I 


where : 




( A. 


- A w> 


n 



is obtained from the correlation curve of " R g W " 


vs Q in Figure 148 for 


Q=R COT 
e LER 




2 

cos 



is the leading- edge radius Reynolds number based on 
LER streamwise values of velocity and radius determined 
at the spanwise location of the mean geometric chord 
of the basic wing. 


and 


A F ‘ A W 


is obtained from Figure 149. 


The upper bound of Region 1 is defined as , which is a func- 
tion of planform geometry and Reynolds number. The planform 
effect for 26.25 million per meter (8.0 million per foot) unit 
Reynolds number on the small-scale models (the "reference 
values") is presented in Figure 150. 


The methodology for determining the other angle-of- attack 
boundaries and the Reynolds number effect on all the angle-of- 
attack boundaries requires that "reference values" of OL^ > Ot 3 , 
# 4 , and ^5 corresponding to 26.25 million per meter (8.0 
million per foot) unit Reynolds number on the small-scale 
models be computed first. Then the Reynolds number corrections 
can be applied to the "reference values". 
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Thus , 


Otr, is obtained from Figure 150. 
ref 


a. = a 9 + 

ref ref 

Aa (3-2) 

from Figure 

a + 

4 ref 3 ref 

A «(4-3) 

from Figure 

a 5 = a 4 + 

ref ref 

Aa (5-4) 

from Figure 


151. 

152. 

153. 


The incremental effect due to Reynolds number requires use 
of two figures (154 and 155) in which an incremental effect of 
Reynolds number is first evaluated for the basic wing planform 
and the effect of irregular planform is applied as a factor. 
Note that in Figure 154 the Reynolds number parameter Q ^ is 
based on the leading-edge radius at the wing tip. 


The effect of Reynolds number is applied to the (X boundary 
values obtained for 26.25 million per meter (8.0 million per 
foot) and model scale, i.e., 


a 2 a 2 f + 
ref 


- _ 


Sa 2 ( A F ) 

A ( Set 2 cos A w ) 


^^2 ^ 


1 

cos X 


W 


where 


A ( 8oc 2 cos a w ) = ( Sa 2 cos a w ) q -(8a 2 cosA w ) & T 


ref 


Thus, for a particular design, Q must be evaluated at model 

ref . . 

scale and at a unit Reynolds number of 26.25 million per meter 
(8.0 million 'per foot), and Q T is evaluated at full scale or 
other desired conditions . 

The other angle-of-attack boundaries (X and OC ^ are 

then calculated in the same way. 
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1 

cos A ^ 

1 _ 

cos A 

w 

1 

cos A w 

The variation of "R(OL)" with « is then constructed by 
sequentially applying the slope functions shown in Figures 156 
through 159 in each appropriate angle-of -attack region as 
follows : 

”r (a)'\ = M R 1 " ; o < a < 

"R (a )" 2 = "R(a )" 1 + (f§^ (a - a 2 ) ; a 2 < a<cc 3 

"R(a )" 3 = "RCa)"! + (f) 2 tas-cQ + (§1 

a 3<a<«4 

"R( a >" 4 = ’TKa)"! + (h|\ (« 3 -«2 ) + (Hi (o: 4-«3 ) 

+ W 4 (a ‘ a4) ; a ^a<«5 

"R(a )" 5 - "RCa)*'! + ( 5 ^ + (H) 3 (a 4* a 3 ) 

+ (H) 4 (a 5- a 4> + (H) 5 (C * - «5> 

; a 5 <a < 26 degrees 
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In practice, one computes corresponding values of Cl from the 
lift prediction method and "RCct)" from the above listed equa- 
tions for a series of angles of attack and solves the basic drag- 
due -to -lift equation. 
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Pitching-Moment Characteristics Predictions 


The prediction methodology with respect to pitching moments 
provides three types of information. First, the aerodynamic- 
center location at low angles of attack can be calculated in 
terms of percent root chord or in terms of the mean geometric 
chord from relatively simple curve fits of data obtained at unit 
Reynolds number of 26.25 million per meter (8.0 million per 
foot) . 


Second, the pitch-up /pitch- down tendencies of the irregular 
planforms at high angles of attack can be assessed on the basis 
of the a.c. location at low angles of attack calculated in terms 
of percent effective mean geometric chord. 

Third, the complete variation of pitching moment (relative 
to 25 percent of the effective mean geometric chord) with angle 
of attack can be constructed. 

Aerodynamic Center In Primary Slope Region (0.1 < < .5) 

The prediction method for aerodynamic center in terms of 
percent root chord is based on curve fits of the experimental 
data for unit Reynolds number of 26.25 million per meter (8.0 
million per foot) . The plot used (Figure 160) presents the 
variations of a.c. location with basic wing leading-edge sweep 
and parametric values of fillet leading-edge sweep. 

A cross-plot of the a.c. location variation with fillet 
leading-edge sweep for a basic wing leading- edge sweep of 45 
degrees was fitted with a 4th degree Legendre polynomial. The 
effects of wing sweep -are accounted for by different linear 
functions in four regions represented by combinations of wing 
and fillet sweep. 

The baseline polynomial for variation of a.c. with 

for = 450 . 
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• (a. c . 'X.cpg = 30.10775858 + 0.1524242511 A p 

+ 0.007543290783 A 2 . 0.0001686738033 A ,, 3 

r r 

+ 0.000001609816225 A p 4 

• To account for other A__ 

W 

• •A < 45°; A p from 25° to 65° 

a.c. 7.^ = (a.c.%C 1 ) B + 5 -^ ( A y - 45°) 

w 

= (a.c.%C 1 ) B 4- (A w - 45°) ( .5700-0.004307692308 A ) 
••A < 45°; A from 65° to 80° 

W — r 

a.c.%C l = (a.c.%C 1 ) B +(A w -45°) [. 2900- . 007733333333 (A F -65°)J 
••A > 45°; A f rom 45 0 to 65° 

a.c.%C 1 =(a.c.%C 1 ) B +(A w -45°) [.2380 + .01035 (A F -45°)J 

• •A., > 45° ; A from 65° to 80° 

a.c.%C 1 =(a.c.%C 1 ) B + (A w -^5°) [.4450 - 0.017 (A p -65°)] 

A comparison of the predicted and test values of aerodyna- 
mic center location is presented in Figure 160. A similar 
approach was used to establish a prediction of the aerodynamic 
center in terms of the percent of mean geometric chord of the 
irregular planform. 

The baseline polynomial for variation of a.c. 7> c ef f with 

A f for A w = 45° is: 

» 

(a.c. % c eff) Ba . = 28.61247819 - .9174258792 A p 

A II =45° +.04500815340 A^ . .000803241088 A 3 

W r . r 

+.000005350951339 A^ 
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To account for other A 

. w 

• A w < 45° ; A p FROM 25° to 60° 

a.c. X c eff = (a. c . 7. c eff > BASIC +( A w *45°) ( . 275- . 00041666666A p ) 

. A„ < 45° ; 

a.c. % c eff 

• A w > 45° ; 

a.c. X c eff = (a.c. 7. c eff ) BASIC +( A w -45°) 

£ 125 + .0078181818 (A p - 45)1 

• A w >45° ; A FROM 72.5° to 80° 
a.c. X c eff = (a.c.% c eff) BASI „+( A w -45°) 

.340 - . 0286666666C A p - 72 . 5)1 

Predicted and test values are compared in Figure 161. 

Pitch-Up/Pitch- Down Tendencies 

Two forms of design guides are presented. The first, 

Figure 162, is a simple matrix presentation showing the tenden- 
cies for each combination of wing sweep and fillet sweep. While 
this guide is useful for interpreting the SHIPS data base, it 
is not necessarily appropriate for analyzing data for other 
configurations. The second guide, Figure 163, makes use of the 
aerodynamic-center location in the primary slope region in terms 
of percent mean geometric chord. For the SHIPS data base it can 
be used with Figure 161 to assess any pitch-up/pitch-down ten- 
dencies. The chart is in the form of a.c. as a percentage of 
effective mean geometric chord rather than percent root chord, 
because the MGC is a more descriptive property of the planform 
than is the root chord. 


A FROM 60° to 80° 

= (a.c.% c eff) BASIC +( A w - 45°) |\ 250- . 0024(A p -60)J 

A„ FROM 45° to 72.5° 

F 
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Pitching 7 Moment Variation With Angle of Attack 

The prediction method for obtaining the variations of 
pitching moment with angle of attack is directly related to the 
drag-due-to-lif t prediction methodology. It was noted that 
breaks or kinks in the pitching -moment curves and axial-force 
curves occurred at or close to the angle-of-attack boundaries 
established in the drag-due-to-lif t analysis. A more extensive 
study showed that two additional boundaries were required, one 
at low angle of attack (Or-.), Figure 164, and one at high angle 
of attack (tt^) , Figure 165. Both of these boundaries have 
corresponding points in the drag- due- to- lift analyses, but they 
were not considered significant because the drag effects are 
very small. For computational purposes, the data presented in 
Figures 164 and 165 were simplified or enhanced to produce a 
useable set of calculation curves. 


As in the drag-due-to-lif t methodology, the pitching-moment 
curves are assumed to be linear functions of angle of attack in 
each region. Except for some irregular planforms which have 
low values of wing sweep and fillet sweep, the linear assumption 
provides a quite accurate representation of the actual curves . 


The methodology is as follows : 
C 
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= c + c o n + c (a, -a..) 
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where : 


OC and Ot ^ are obtained from Figures 164 and 165 , and the 


values of C 


m 


in 




nt 


m 


a. 


a=o ‘“a 0 “a^ n bt 2 ’ “tar,’ “a 4 

are obtained from Figures 166 through 172. 


C , and 

a. 


The variation of pitching moment at zero angle of attack 
presented in Figure 167 represents smoothed average values for 
all outboard panel sweeps . The maximum deviation from the curve 
was + 0.004 at A ^ = 65, which is well within the desired accuracy. 


Note that for the NACA 0008 airfoils the values of C 


m Q^ 


and C are identical, but provision has been made in the 


m 


a 


methodology to account for the fact that for thicker airfoils 

C and C may have different values, 

nu, nu, J 


Since the intent of the pitching- moment calculation is pri- 
marily to define the stability characteristics of the configura- 
tion, the basic prediction method does not account for Reynolds 
number effects on either the angle -of -attack boundaries or on 
the pitching-moment slopes because the lift prediction does not 
include such effects. The effect of Reynolds number on the 
pitching-moment slopes is small for the NACA 0008 airfoils , but 
some effect was noted for the thicker airfoils from the basic 
data plots. These effects were to be accounted for in the final 
methodology . 
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Figure 146. Lifting Surface Correlation Factor for Subsonic 
Minimum Drag 
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Figure 149. Calculation Chart for Incremental Effect of Fillet 
Sweep on Suction Ratio in Region 1 
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Figure 151. Incremental Angle of Attack for Region 2 
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Figure 162. - Pitch Up/Pitch Down Tendencies as 
Functions of Wing and Fillet Sweep 
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Figure 169. Pitching -Moment Coefficient Derivative for 
Region 3 
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Figure 170. Pitching Moment Coefficient Derivative for 
Region 4 

















REVISIONS TO "BASIC" METHODS 


The basic prediction methods presented in the previous section 
partially achieved the objectives of the investigation. Ongoing 
NASA configuration definition studies of an advanced aerospace 
vehicle indicated the need for data on configurations with other 
airfoils than the NACA 0008 airfoil. The very limited data ob- 
tained with three configurations having NACA 0012 airfoils indi- 
cated that Reynolds number effects might be significant to much 
larger Reynolds numbers than had been obtained in test ARC 12-086. 
In addition, limitations on the angles of attack obtained in the 
first test had left some questions as to the pitching -moment be- 
havior for some configurations. Thus, it was apparent even before 
the "basic" prediction method task was completed that additional 
testing would be required to provide data to resolve these is- 
sues . 


Additional Testing and Analysis 


The Vehicle Analysis Branch of the NASA Langley Research Cen- 
ter planned the additional tests of a limited series of plan- 
forms (described earlier in the Test Data Base section of this 
report) to define the effects of different airfoil sections and 
higher Reynolds numbers, and eventually data were obtained at 
higher angles of attack. The tests were accomplished in the NASA 
Langley Research Center Low Turbulence Pressure tunnel with the 
small models and the NASA Ames Research Center 12- foot tunnel with 
the large (twice-size) models. 

Analysis of the additional data was performed using the same 
approaches that were used for the "basic" prediction method de- 
velopment. This required that the reference quantities of wing 
area, moment reference length, and moment reference point be on 
the same basis, whereas the test data for the additional tests had 
been reduced originally using different quantities. The neces- 
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sary revisions to the data were made by altering a few equations 
in the analysis procedures which had been programmed for solution 
on a Hewlett Packard 9820 desktop computer. 

One of the first steps in the analysis was to make compari- 
sons of data obtained in the different facilities to determine 
if significant differences occurred. An example comparison of 
the variations of pitching moment with angle of attack is pre- 
sented in Figure 173. It is apparent that the data for the 
65A012 airfoil configuration obtained in LTPT test 262 is signi- 
ficantly different from that for all of the other configurations. 
In contrast, data from test ARC 12-257 shows only small difference 
between the pitching-moment curves for 0012 and 65A012 airfoil 
configurations. It later became apparent that erroneous moment 
transfer distances had been used in reducing data for four con- 
figurations in test LTPT 262, including the 70 and 75 degree fil- 
let configuration with 65A012 and 65^412 airfoils. It became 
necessary later to adjust the moment- transfer distances for 
these configurations in order to complete the pit chin g- moment pre 
dictions. The adjustments were made such that the slopes at low 
angle of attack agreed with the data from ARC test 12-257. 

An overall survey of Reynolds number effects indicated that 
above a unit Reynolds number of 26.25 million per meter the ef- 
fects are small in comparison to the effects noted between 13.13 
million per meter and 26.25 million per meter for the small mo- 
dels. The configurations with 12 percent thick airfoils show 
the most sensitivity to Reynolds number. It is apparent that the 
major effects occur as a result of relatively low leading-edge 
sweep of the outboard wing panel or of the fillet-outboard panel 
combination. 

A brief look at the suction- ratio data showed that, in gen- 
eral, the values in region 1 are higher than those measured in 
ARC test 12-086. The reason for this fact is not immediately 
obvious. It was therefore decided to concentrate first on the 
lift correlation to see if significant differences occurred bet- 
ween the tests in the different facilities and with the differ- 
ent size models and to look at the effects of airfoil thickness 
ratio, thickness distribution, and camber. 
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Some example comarisons are presented in Figures 174 
through 179. Figures 174 and 175 illustrate the magnitudes of 
the differences in the variation of the lift- correlation 
parameter with angle of attack for three of the 4 tests in which 
the same configurations were run. Also shown are the predicted 
values obtained from the overall correlation procedure. It is 
apparent that measurements made in the Low Turbulence Pressure 
Tunnel on the small models and in the ARC 12- ft tunnel on the 
large models produced higher lift coefficients at a given angle 
of attack than did the original measurements in the ARC 12- foot 
tunnel on the small models. The differences are of the order of 
4 to 7 percent in the regions where the curves are reasonably 
smooth. It is not at all obvious that any one set of data is 
more correct than the others. 

Figure 176 and 177 present the effects of airfoil- thickness 
ratio variations as measured in the LTPT and ARC 12-ft facilities 
on the small and large models of Wing III with the 80 degree fil- 
let. Although there are differences in details of the curves, the 
general trends of the effects are similar. The same comment ap- 
plies to the data for the effects of airfoil-thickness ratios 
presented in Figures 178 and 179 for Wing III with no fillet. 

Figures 180, 181 and 182 show the effects of airfoil thick- 
ness distribution and camber for the three planforms having 
Wing III with 80, 75 and 70 degree fillets, respectively. The 
trends with fillet sweep are consistent. The camber effect de- 
creases with increasing fillet sweep. The large model data pro- 
duces similar trends. 

Figures 183, 184 and 185 show the effect of Reynolds number 
for Wing III with 70 degree fillet and the three 12 -percent thick 
airfoil sections. Note that the effects are primarily in the 
stall progression region and, therefore, would affect the stall 
progression factor and the incremental lift charts. 

One approach considered for modifying or revising the pre- 
diction method's was to obtain incremental lift curve parameter 
values from the data obtained on the NACA 0008 airfoil configura- 
tions at unit Reynolds number of 26.25 million per meter (8.0 
million per foot). Figures 186-190 present the effect of fillet 
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sweep on the lift curve parameter increments as functions of 
angle of attack for the different airfoil sections as tested in 
the LTPT facility. For each airfoil, a reasonably consistent 
trend exists for the effects of fillet sweep. However, when the 
effects of Reynolds number are included as in Figure 191 the pic- 
ture becomes confused. 

The approach of trying to use the NACA 0008 airfoil data as 
a baseline for the airfoil effects and obtain increments for the 
airfoil effects and other increments for the Reynolds number ef- 
fects was not successful. There were too many conflicting 
trends in the data to evolve simple correlations. As a conse- 
quence, it was decided to use a more direct approach. 

The plots of Reynolds number effects on lift, suction ratio, 
and pitching moment were reviewed one more time and it was noted 
that the key results could be represented by defining the charac- 
teristics for three values of unit Reynolds number: (1) the 

minimum value tested (13.13 million per meter or 4.0 million per 
foot), (2) a "standard" value corresponding to the maximum value 
tested in the original planform series (26.25 million per meter 
or 8.0 million per foot), and (3) the maximum value tested in the 
additional test series (Usually 45.93 million per meter or 14.0 
million per foot). 

A prediction using the minimum Reynolds number data could 
be used to help define loads for wind-tunnel models to be tested 
in facilities having only relatively low unit Reynolds number 
capabilities . 

The use of the "standard" Reynolds number data would allow 
a tie-in between the airfoil effects and the planform effects 
obtained from the "basic" prediction method. For practical pur- 
poses the maximum unit Reynolds number data might be considered 
representative of what could be expected for full-scale flight. 

This approach was taken and a set of revised prediction 
methods were developed. The revised methods use slightly dif- 
ferent chart formats and equations for two reasons : 
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(1) The maximum angle of attack was extended from 26 de- 
grees to as high as 32 degrees. 

(2) The limited series of planforms tested consisted 
only of Basic' Wing III with various fillets. Thus 
the effect of outer panel sweep could not be con- 
sidered in the revised methods. 


Revised Lift Prediction Method 


The revised prediction method for the lift characteristics 
uses the same correlation parameter as the basic prediction meth- 
od, but in a slightly different way. The revised lift prediction 
equation is: r , . 1 


C 


L 





std 



RN 


( )J 


The terms in this equation are the same as for the basic prediction 
method with the following exceptions: 


(1) The terms in brackets are defined by sets of charts 
appropriate to a specific airfoil section. 


(2) The first term in brackets is obtained from a correla- 


tion chart which has the lift parameter 


against the parameter 
angle of attack. 



'std 


/StanAp for parametric values of 


The values obtained from this term represent data ob- 
tained at the Standard unit Reynolds number of 26.25 
million per meter (8.0 million per foot). 


(3) The second term in brackets represents the incremental 
effect of Reynolds number and is obtained from charts 
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which show the incremental value 



for 


either RN . or RN 


„ mm , max 

angles of attack. 


in carpet plot form for various 


Table 4 presents a key to the Figures which contain 
the sets of charts for each airfoil. 


TABLE 4 KEY TO REVISED LIFT PREDICTION CHARTS 


AIRFOIL 

SECTION 



std 


C V 

L tt b /rn 


min 


j a 



max 


NACA 0004 

Fig 192 

N. A. 

N. A. 

NACA 0008 

Fig 193 

Fig 194 

Fig 195 

NACA 0012 

Fig 196 

Fig 197 

Fig 198 

NACA 0015 

Fig 199 

Fig 200 

Fig 201 

NACA 65A012 

Fig 202 

Fig 203 

Fig 204 

NACA 65 x 412 

Fig 205 

Fig 206 

Fig 207 

THIN HEXAGON 

Fig 208 

N. A. 

N. A. 

WITH SHARP L.E. 



. 

Figures 192, 

193, 196, and 

199 present 

the lift corre 


curves for the NACA 0004, 0008, 0012 and 0015 airfoils, respect- 
ively. The NACA 65A012, 65^412, and thin hexagon airfoils were 
only tested for a few conf igurations , Figures 202 and 205 show 
comparisons of the lift correlation curves of the 65A012 and 
65,412 airfoils tested with the corresponding curves for the 
NACA 0012 airfoil, and Figure 208 shows comparison of the lift 
correlation curves of the 4 percent thicK hexagon airfoils 
with corresponding curves for the NACA 0004 airfoils. If for 
some reason it were necessary to estimate the lift curve of a 
configuration having the 65A012 or 65^412 airfoils with a fillet 
sweep less than 70 degrees,, the appropriate correlation curves 
could be extended to lower values of/3tanA F using the NACA 0012 
curves as a guide. 

With respect to the incremental values due to Reynolds 
number, it was found that for the 4-percent thick airfoils 
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the effects of Reynolds number were not significant. 

The prediction method consists of calculating values for 
several values of angle of attack using the appropriate set of 
charts for the desired configuration and airfoil section. The 
charts are set up to produce lift curves for any thickness ratio 
between .04 and .15. 

In order to provide a tie-in to the "basic” prediction 
method, data for the Wing III series of planforms with NACA 
0008 airfoils from test ARC 12-086 were plotted in the format 
used for the revised prediction method. The resulting chart 
is presented in Figure 209. 


Figure 209 can be of value in two ways. First, it can be 
used to establish the uncertainties in lift predictions for 
the family of irregular planforms having Wing III as the basic 
planform, by comparing results obtained using Figure 209 and 
Figure 193. Second, it provides a means of applying the incre- 
mental effects of airfoil thickness ratio, airfoil shape, and 
Reynolds number determined from the revised prediction method 
to a broader range of planforms. It is likely. that reasonable 
estimates of the lift curves could be made for irregular plan- 
forms having basic planforms in the range from Wing II to 
Wing IV, by applying the increments to results obtained from 
the "basic" prediction method for the desired planform. 
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Revised Drag- Predict ion Method 


Revisions to the drag, prediction methodology affected only 
the drag-due-to-lift part of the method. Revisions were re- 
quired for three reasons: 

(1) The plateau values of suction ratio for the basic 
wings of different thickness ratios, as tested in 
the LTPT facility, were significant ly higher and 
formed a different correlation curve as a function 
of the Q parameter from the curve developed for the 
basic prediction method. 

(2) The slopes of the suction-ratio curves in the var- 
ious regions varied significantly with Reynolds num- 
ber for the thicker airfoils for some of the plan- 
forms . 

(3) The approach used to predict the effects of Reynolds 
number on lift was not compatible with the method 
used to account for Reynolds number effects on drag 
due to lift in the original formulation used in the 
basic prediction method. 

The revised method still makes use of the suction- ratio 
concept and the representation of the variation of suction ratio 
with angle of attack by a series of linear segments as illu- 
strated in Figure 210. The angle-of -attack boundaries defining 
the segments also apply to the pitching-moment prediction and 
in fact some iteration between values appropriate to the suc- 
tion-ratio curves and those appropriate for the pitching- moment 
curves was required to produce compatible results. 
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As many as 7 linear segments were required to produce rea- 
sonable fit's to the experimental suction ratio and/or pitch- 
ing moment data because the angle-of-attack range for the pre- 
diction was extended from 26 to 32 degrees. Some of the con- 
figurations did not require 7 segments. In those cases, ficti- 
tious boundaries were defined to allow a single generalized 
equation to represent all the configurations. 

The basic drag-due -to-lift equation is the same as used 
for the "basic" prediction method. 


C D t C tana- "R(a) M 
L L 

but for the revised method 


C^tan a - 


L 
7T A 


"R(a) M = R- 


; o <a<a 


Rf + 


d"R"\ (a-a->) ;a 9 <a<a 

da L 


r i + 


d"R' 

d a 


(a, -a,) + / d"R" \ (a-oL) 
l d “ 3 


; a 3 < oc < a 4 


= \ + (0 (« 3 -«2 ) + «V<V 


+ /d M R M \ 


; a 4 <a<a 5 
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( d n R M 

V da , 

) < a 3 -<V + | 
2 

/d"R"\ 

l dot 4 

«V°3> 

/ d"R"' 

\ d a ; 

) (a 5 -a 4 ) + 

/ d"R”\ 

v da y 

(a-a 5 ) 


4 

\ / 

5 



: oc 5 ^a<a 6 


d"R' 
d a 


^ «X 3 - rtj) + «v a 3> 


r d"R' 
d a 


j (a 5 -a 4 ) + ^ d M R M j f<* 6 -a 5 ) 


d n R" \ (a-<x ) 

doTj a 6 
'6 


; a 6 <a<a7 


/ d"R" \ 

v da A 


(a 3 -a 2 ) + «V« 3 ) 


3 


d"R" \ «X 7 -a 6 ) + 


d a 



(a-oc 7 ) 


; a 7 <a< 32 

The values of the CL boundaries , the plateau value of the suc- 
tion ratio, R, , and the values of the slopes of the suction ra- 

/d ,, R n \ 

tio curves in each region, | J were determined for each con- 
figuration for each of the three Reynolds number conditions 
used in the revised lift prediction method. The CL boundaries 
are tabulated in Tables 5 through 8, and the plateau values of 
suction ratio and the slopes of the various segments of the 
suction ratio variations with angle of attack are tabulated in 
Tables 9 through 12. 


To use the prediction method, the calculations are made for 
a series of angles of attack using compatible values for the 
various terms, i.e., the planform, airfoil section, and Reynolds 
number condition must be the same for both lift and the drag- due 
to-lift conditions. If the user wishes to predict full-scale 
flight characteristics, the lift and drag-due-to-lift terms for 
R^ conditions should be used. The minimum-drag term should 
max 

be calculated for the actual flight Reynolds number condition. 


Revised Pitching- Moment Prediction Method 


The only- revisions to the pitching-moment prediction meth- 
od are those required to extend the angle -of-at tack range from 
26 to 32 degrees and to make the method consistent with the re- 
vised lift and drag-due-to-lift predictions. 
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The. equation for the variation of pitching moment with 
angle of attack is: 


= c + c (a -o) ; o<a<a, 

rr> 1 


m . 25c m a=o “a 


o 


= c 


„ + V. (a r 0) + C m_. 


a=o l a- 


ct. 


= c + c (a.-o) + c (a 0 -OL) 

m xru ( 1 m^y 2 i 


a=o «o 


+ c 


(a ^2> 


«2 


} at 2 <oc<a 3 


C m + C m (a r 0) + C m 

a=o “o “1 


+C rn «W + C m„ \ ; tt 3* 

2 a 3 


c + c (<x.- o) + c (a 0 -a,) 
m nu, 1 m_, 2 

a=o a o °i 


+ c (aL-aJ + c (a A -a,) 

nu, 3 2 m 4 3 

W 0 3 


+ c (a-a, ) 
m a f 4 

4 


. a 4 <a<a^ 
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(a 2 -ap 


c - c + c (a,-o> + c 

m . 25c "b'O 1 


+ V <a 3- a 2> + C m (a 4 ' a 3 ) 
z 3 


+ c (a, -a ) + c («-<x) ; a 

HU, J -4 HU, ^ j 

4 5 


- c + c (a,-o) + c (a, -a.) 
"to-o a 0 X 


+ V (a 3 '«2> + (a 4' Q 3 ) 


a, 


+ c («,-«,) + c (a, -at) 

X X 


+ c 


m, 


a. 






= c 


m + V «v°> + V *W> 


a=o 


°b 


a. 


+ C m ( °W + «V<V 


a 


a , 


+ c (a, -a,) + c ( Ptr-CL ) 

nu 5 4 nu, ° 


'+ c (a,-<x) + c 

X 1 «7 


<a<a 6 


a<a 7 


(a-a ? ) ; a 7 <o:<32 



The angle -of- attack boundaries used for this equation are iden- 
tical to those used for the suction-ratio calculation which 
were presented in Tables 5 through 8. 

The values of C and the slopes of the various seg- 

m a=o 

ments of pitching-moment variations with angle of attack are 
presented in Tables 13 through 16 for the various combinations 
of* planforms, airfoil section, and Reynolds number condition. 
Obviously, the pitching-moment prediction requires use of the 
same planform, airfoil section, and Reynolds number conditions 
as are used for the lift and drag-due- to-lift calculations. 

The "revised" prediction methods are strictly applicable 
only to the Wing III series of irregular planforms with fillet 
sweeps between 45 degrees and 80 degrees, whereas the "basic" 
prediction methods are applicable to a wide range of irregular 
planforms having NACA 0008 airfoils but only at the standard 
Reynolds number condition. 

It is considered feasible to apply judicious combinations 
of the two sets of prediction methods to obtain the lift, drag, 
'and pitching-moment variations with angle of attack for configu- 
rations which have planforms which are not of the Wing III 
series but have airfoils other than the NACA 0008 airfoil sec- 
tion. 
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TABLE 5 


Angle-of-Attack Boundaries for Drag Due to Lift and 
Pitching-Moment Segments - NACA OOXX Airfoils, Rm 

min 


CONFIGURATION 

A A 

oc 

1 

“2 

Of 

3 

CL 

4 

a 

5 

Of 

6 

“7 

F 


AIRFOIL 








45 

45 

0004 

2.4 

2.4 

7,2 

13.4 

22.7 

26.6 

30.5 



0008 

4.4 

4.4 

7.5 

15.4 

20.6 

27.8 

31.2 



0012 

6.0 

6.0 

10.8 

10.9 

17.5 

22.0 

26.0 



0015 

7.0 

7.0 

12.0 

17.7 

19.8 

25.8 

30.0 

60 

45 

0004 

3.0 

3.0 

6.35 

15.6 

19. 3 

23.5 

30.6 

65 

45 

0004 

3.0 

3.0 

7.0 

15.6 

19.3 

23.5 

30.6 



0008 

5.6 

5.6 

10.2 

14.2 

17.0 

22.4 

27.0 



0012 

7.4 

7.4 

13.5 

15.6 

18.0 

22.6 

29.0 



0015 

6.0 

6.0 

12.2 

13.3 

20.0 

24.0 

28.0 

70 

45 

0004 

3.6 

3.6 

6.8 

11.5 

15.8 

18.8 

28.3 



0008 

6.6 

6.6 

10.0 

14.2 

20.2 , 

24.6 

28.5 



0012 

7.0 

7.0 

13.5 

18.0 

22.5 

26.6 

30.5 



0015 

5.6 

5.6 

14.0 

17.2 

20. 6 

27.2 

30.0 

75 

45 

0004 

3.0 

3.0 

7.1 

12.1 

18.4 

25.6 

29.0 



0008 

4.5 

4.5 

7.1 

11. 3 

17.1 

20.8 

26.4 



0012 

6.0 

6.0 

8.6 

13.2 

17.5 

21.1 

25.0 



0015 

4.4 

4.4 

8.0 

13. 8 

17. 9 

20.8 

30.0 

80 

45 

0004 

2.6 

2.6 

7.95 

14.5 

17.3 

25.5 

29.3 



0008 

4.8 

4.8 

7.5 

12.5 

18.1 

23.9 

29.3 



0012 

6.4 

4.8 

11.1 

15.7 

18. 7 

25.6 

30.0 



0015 

7.0 

7.0 

12.0 

15.7 

20.0 

23.0 

28.0 



TABLE 6 


Angle-of- Attack Boundaries for Drag Due to Lift and 
Pitching-Moment Segments - NACA OOXX Airfoils, % st< j 


CONFIGURATION 

a f a w airfoil 

*1 

« 2 

tt 3 

*4 

“5 

«6 

*7 

45 45 

0004 

2.4 

2.4 

8.4 

13.8 

22. 7 

26.6 

30.5 


0008 

4.4 

4.4 

9.6 

16.9 

23.0 

27.8 

31.2 


0012 

6.0 

6.0 

10.8 

14.4 

21.0 

25.6 

31. 1 


0015 

7.0 

7.0 

12. 0 

19.1 

21.0 

24.2 

28.0 

60 45 

0004 

3-0 

3. 0 

6. 6 

11. 7 

14. 8 

19.2 

31.2 

65 45 

0004 

3. 0 

3.0 

7.0 

15.6 

19.3 

23.5 

30.6 


0008 

5.6 

5.6 

11.2 

16.0 

19.6 

24.4 

28.0 


0012 

7.4 

7.4 

11. 7 

16.0 

20.0 

25.2 

32.0 


0015 

6.0 

6.0 

12.2 

16.0 

19.8 

24.0 

26.4 

70 45 

0004 

3.6 

3.6 

7.2 

12.7 

15.8 

18.8 

28.3 


0008 

6.6 

6. 6 

11. 0 

17. 0 

21. 8 

24. 6 

29.6 


0012 

7.0 

7.0 

13.9 

17.5 

21.5 

25.6 

30.0 


0015 

5.6 

5.6 

14.0 

18.0 

21.0 

25.6 

30.4 

75 45 

0004 

3.0 

3.0 

9.4 

12. 1 

18.4 

25.6 

29.0 


0008 

4.5 

4. 5 

8. 1 

13.8 

19. 8 

23. 7 

28.0 


0012 

6.0 

6.0 

8. 1 

14. 1 

17.8 

22.0 

28.0 


0015 

4.4 

4.4 

8.0 

13.8 

17.9 

20.8 

30.4 

80 45 

0004 

2.6 

2.6 

8.8 

14.5 

17.3 

25.5 

29.3 


0008 

4. 8 

4. 8 

9.5 

12.6 

17.4 

23.9 

28.8 


0012 

6.4 

6 . 4 

10.8 

12.2 

17. 3 

24.0 

31.2 


0015 

7.0 

7.0 

10. 8 

12.0 

20.0 

27.2 

33.0 



TABLE 7 


Angle-of- Attack Boundaries for Drag Due to Lift and 
Pitching- Moment Segments - NACA OOXX Airfoils, % max 


CONFIGURATION 

A A 

*1 

«2 

a 3 

«4 

a 5 

tt 6 

tt 7 

a F 

A W 

AIRFOIL 








45 

45 

0004 

2.4 

2.4 

9.2 

13.8 

22.7 

26.6 

30.5 



0008 

4.4 

4.4 

12.2 

16.0 

18.4 

23.0 

27.8 



0012 

6.0 

6.0 

12.6 

20.6 

22.2 

23.5 

27.5 



0015 

7.0 

7.0 

12.0 

18.2 

22.2 

23.5 

29.5 

60 

45 

0004 

3.0 

3.0 

7.6 

11.7 

14.8 

19.2 

31.2 

65 

45 

0004 

3.0 

3.0 

7.0 

15.6 

19. 3 

23.5 

30.6 



0008 

5.6 

5.6 

11.15 

16.0 

19.6 

24.4 

28.0 



0012 

7.4 

7.4 

13.0 

17.9 

20.8 

28.6 

32.0 



0015 

6.0 

6.0 

12.2 

14.2 

18.0 

20.3 

25.2 

70 

45 

0004 

3.6 

3.6 

7.65 

12.7 

15.8 

18.8 

28.3 



0008 

6.6 

6.6 

11.65 

16.7 

21.0 

23.6 

27.0 



0012 

7.0 

7.0 

14.7 

17.7 

21. 8 

25.3 

30.0 



0015 

5.6 

5.6 

14.8 

17.6 

19.3 

27.0 

31.0 

75 

45 

0004 

3.0 

3.0 

9.8 

12.1 

18.4 

25.6 

29.0 



0008 

4.5 

4.5 

8.0 

12.45 

19.8 

24.0 

29.0 



0012 

6.0 

6.0 

8. 1 

16.6 

20.5 

24.5 

28.6 



0015 

4.4 

4.4 

12.0 

16.9 

20.1 

25.6 

29.6 

80 

45 

0004 

2.6 

2.6 

9.1 

14.5 

17. 3 

25.5 

29.3 



0008 

4.8 

4. 8 

9.5 

13.1 

16.9 

22.5 

29.1 



0012 

6.4 

6.4 

9.8 

16.0 

18.9 

25.0 

28.0 



0015 

7.0 

7.0 

10.8 

12.0 

20.0 

24.4 

29.0 



TABLE 8 


Angle -of- At tack Boundaries for Drag Due to Lift an* 7 
Pitching-Moment Segments, NACA 65A012 and 65x412 


CONFIGURATION 
A F A W AIRFOIL 

RN 

®1 

a 2 

«3 

«4 

ol 

j 

«6 

«7 

70 

45 

65A012 

Min. 

7.4 

7.4 

9.6 

15.7 

20.3 

24. 8 

27.5 

75 

45 

6 5 A0 12 

Min. 

5.4 

5.4 

9.1 

11. 1 

16.9 

22.0 

26.1 

80 

45 

6 5 A0 12 

Min. 

7. 1 

7.1 

12.2 

15.4 

20.2 

26.0 

31.4 

70 

45 

6 5 A0 12 

Std. 

7.4 

7.4 

10.2 

14.8 

20.4 

24.6 

30.8 

75 

45 

6 5 A0 12 

Std. 

4 . 6 

4.6 

9.6 

17.0 

20.8 

24.6 

27.8 

80 

45 

65A012 

Std. 

5.8 

5.8 

10. 8 

15. 1 

20.0 

27.2 

31.0 

70 

45 

65A012 

Max. 

7.4 

7.4 

10.2 

15.8 

20. 1 

25.0 

28.6 

75 

45 

65A012 

Max. 

4. 8 

4. 8 

10.0 

16.2 

18.4 

21.6 

28.2 

80 

45 

6 5 A0 12 

Max. 

6.0 

6.0 

12.0 

15. 8 

21.2 

27.7 

32.0 

70 

45 

65]412 

Min . 

7.4 

7.4 

10. 7 

13.8 

17.7 

22.8 

27.2 

75 

45 

65i412 

Min. 

6.9 

6.9 

14.4 

18.2 

20.7 

27. 3 

32.0 

80 

45 

65 i412 

Min. 

5.6 

5.6 

13.0 

18.9 

23.5 

28.0 

30.5 

70 

45 

65 X 412 

Std. 

7.6 

7.6 

10.3 

14.1 

19.1 

25.6 

32.4 

75 

45 

65x412 

Std. 

6.9 

6.9 

15. 3 

17. 8 

22.4 

28.4 

31.0 

80 

45 

65x412 

Std. 

5.5 

5.5 

10.8 

15.4 

20.6 

25.7 

30.1 

70 

45 

65x412 

Max. 

7.4 

7.4 

11.0 

15.4 

19.2 

22.6 

28.6 

75 

45 

65x412 

Max. 

6.9 

6.9 

13.6 

16.2 

19.2 

24.0 

29.0 

80 

45 

65x412 

Max. 

5.6 

5.6 

13.0 

16.0 

22.8 

28.0 

32.0 

45 

45 

00SLE 

Std. 

2.4 

2.4 

7.2 

15. 8 

20. 0 

22.9 

26.9 

75 

45 

OOSLE 

Std. 

2.0 

2.0 

7.0 

13.0 

18.2 

23.4 

31.6 

80 

45 

OOSLE 

Std. 

2.0 

2.0 

7.9 

11.6 

17.2 

23. 8 

28.0 
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TABLE 9 


Plateau Values of Suction Ratio and Slopes of 
Suction-Ratio Segments - NACA OOXX Airfoils, 








dR 

dR 

dR 

dR 

dR 

dR 

CONFIGURATION 


R 1 

da 2 

dC* 3 

da 4 

da 5 

d« 6 

d a 7 

\ 

A w 

AIRFOIL 










45 

45 

0004 

Min. 

.885 

-.0727 

-.0355 

-.0197 

-.0102 

-.0070 

-.0031 



0008 



.967 

- . 0032 

-.0705 

-.0443 

-.0075 

-.0028 

-.0028 



0012 


1 

.980 

-.0040 

-.0040 

-.0874 

-.0488 

0 

0 



0015 



.985 

-.0037 

-.0037 

-.1960 

-.0260 

-.0238 

-.0238 

60 

45 

0004 

Min. 

.838 

-.0893 

-.0414 

-.0290 

-.0158 

-.0088 

-.0036 

65 

45 

0004 

Min. 

.835 

-.1175 

-.0229 

-.0103 

-.0126 

-.0077 

-.0048 



0008 



. 860 

-.0048 

-.0844 

-.0417 

-.0215 

-.0215 

-.0060 



0012 

! 


. 905 

-.0010 

-.1484 

-.0812 

-.0292 

-.0012 

-.0222 



0015 



.940 

-.0039 

-.0039 

-.0470 

-.0817 

-.0018 

-.0018 

70 

45 

0004 

Min . 

.840 

-.1254 

-.0315 

-.0142 

-.0142 

-.0092 

-.0048 



0008 



.950 

-.0498 

-.0846 

-.0102 

-.0245 

-.0168 

-.0087 



0012 



.955 

-.0067 

-.0865 

-.0418 

-.0101 

-.0122 

-.0114 



0015 


t 

.990 

-.0082 

-.0082 

-.0930 

-.0300 

-.0140 

-.0140 

75 

45 

0004 

Min. 

.810 

-.0943 

- .0310 

-.0137 

-.0068 

-.0030 

- .0047 



0008 

I 

i 


.910 

-.0087 

-.0931 

-.0338 

-.0055 

-.0166 

-.0056 



0012 



.955 

-.0022 

-.0186 

-.0638 

-.0475 

-.0092 

-.0182 



0015 

( 


.965 

-.0120 

-.0142 

-.0342 

-.0906 

-.0125 

-.0048 

80 

45 

0004 

Mi 

n. 

. 882 

-.1070 

-.0161 

-.0096 

-.0063 

-.0038 

-.0078 



0008 

1 


.915 

-.0337 

-.0707 

-.0392 

-.0055 

-.0079 

-.0130 



0012 

i 


.912 

-.0327 

-.0327 

-.0593 

-.0016 

-.0016 

-.0016 



0015 

i 


.915 

-.0220 

-.0127 

-.0378 

-.0683 

-.0080 

-.0080 
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TABLE 10 


Plateau Values of Suction Ratio and Slopes of 
Suction-Ratio Segments - NACA 00XX Airfoils, Rm 

™std 


CONFIGURATION 

a f a w airfoil 


R 1 

dR 

da 

2 

dR 

d “ 3 

dR 

tt 4 

dR 

da 

5 

dR 

a 

6 

dR 

da 

7 

A 5 

45 

0004 

Std 

.940 

-.0727 

-.0355 

-.0197 

-.0102 

-.0070 

-.0031 



0008 



.980 

-.0032 

-.0705 

-.0443 

-.0075 

-.0028 

-.0028 



0012 



.990 

-.0040 

-.0040 

-.0721 

-.0512 

-.0055 

-.0061 



0015 



.965 

-.0037 

-.0037 

-.1560 

-.1032 

-.0229 

-.0010 

60 

45 

0004 

Std 

.882 

-.0893 

-.0414 

-.0290 

-.0158 

-.0088 

-.0036 

65 

45 

0004 

Std 

.900 

-.1175 

-.0229 

-.0103 

-.0126 

-.0077 

-.0048 



0008 



.940 

-.0048 

-.0844 

-.0417 

-.0153 

-.0153 

-.0060 



0012 



.925 

-.0010 

-.0010 

-.0932 

-.0451 

-.0163 

-.0310 



0015 

i 


.940 

-.0039 

-.0039 

-.0145 

-.0817 

-.0436 

-.0146 

70 

45 

0004 

Std 

.905 

-.1254 

-.0315 

-.0142 

-.0142 

-.0092 

-.0048 



0008 



.975 

-.0198 

-.0924 

0 

-.0204 

-.0167 

-.0167 



0012 



. 960 

-.0067 

-.0833 

-.0429 

-.0244 

-.0163 

-.0118 



0015 

i 


.980 

-.0082 

-.0118 

-.0609 

-.0431 

-.0150 

-.0218 

75 

45 

0004 

Std 

.930 

-.0943 

-.0310 

-.0137 

-.0008 

-.0030 

-.0047 



0008 



.945 

-.0087 

-.0708 

-.0465 

-.0091 

-.0112 

-.0180 



0012 



.955 

-.0022 

-.0172 

-.0771 

-.0473 

-.0093 

-.0184 



0015 



.965 

-.0120 

-.0142 

-.0342 

-.0906 

-.0125 

-.0048 

80 

45 

0004 

Std 

. 960 

-.1070 

-.0161 

-.0096 

- .0063 

-.0038 

-.0078 



0008 



. 940 

-.0337 

-.0408 

-.0755 

-.0057 

-.0067 

-.0067 



0012 



.912 

-.0262 

-.0262 

-.0304 

-.0427 

-.0027 

-.0027 



0015 

i 


.915 

-.0220 

-.0220 

-.0127 

-.0345 

-.0138 

-.0072 
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TABLE XI 


Plateau Values of Suction Ratio and Slopes of 
Suction-Ratio Segments - NACA OOXX Airfoils, % 

max 


CONFIGURATION 
a f a w AIRFOIL 

*N 

R 1 

dR 

d<X 2 

dR 

d« 3 

dR 

^4 

dR 

d<* 5 

dR 

da, 

6 

dR 

dCX 7 

45 

45 

0004 

Max. 

.970 

-.0727 

-.0355 

-.0197 

-.0102 

-.0070 

-.0031 



0008 



. 980 

-.0032 

-.1048 

-.0623 

-.0443 

-.0075 

-.0028 



0012 



.982 

-.0040 

-.0040 

-.2470 

-.2470 

-.0160 

0 



oo;5 



.985 

-.0037 

-.0037 

-.1385 

-.1385 

-.0090 

-.0010 

60 

45 

0004 

Max. 

.925 

-.0893 

-.0414 

-.0290 

-.0158 

-.0088 

-.0036 

65 

45 

0004 

Max. 

.962 

-.1175 

-.0229 

-.0103 

-.0126 

-.0077 

-.0048 



0008 



.956 

-.0048 

-.0844 

-.0417 

-.0153 

-.0153 

-.0150 



0012 



.925 

-.0010 

-.0055 

-.1144 

-.0225 

-.0150 

-.0150 



0015 

i 


.940 

-.0039 

0.0039 

-.0110 

-.2000 

-.0101 

-.0141 

70 

45 

0004 

Max. 

.945 

-.1254 

-.0315 

-.0142 

-.0142 

-.0092 

-.0048 



0008 



.982 

-.0080 

-.1073 

-.0167 

.0197 

-.0242 

-.0108 



0012 



.965 

-.0048 

-.0829 

-.0587 

-.0233 

-.0150 

-.0218 



0015 


t 

.990 

-.0070 

-.0160 

-.1060 

-.0366 

-.0148 

-.0148 

75 

45 

0004 

Max. 

.972 

-.0943 

-.0310 

-.0137 

-.0068 

-.0030 

-.0047 



0008 



.970 

-.0087 

-.0422 

-.0685 

-.0020 

-.0167 

-.0167 



0012 



.955 

-.0070 

-.0165 

-.0895 

-.0255 

-.0087 

-.0087 



0015 


< 

.980 

-.0120 

-.0120 

-.0178 

-.0925 

-.0214 

-.0022 

80 

45 

0004 

Max. 

1.000 

-.1070 

-.0161 

-.0096 

-.0063 

-.0038 

-.0078 



0008 



.940 

-.0322 

-.0322 

-.0965 

-.0062 

-.0076 

-.0076 



0012 



.912 

-.0262 

-.0175 

-.0990 

-.0320 

-.0076 

-.0076 



0015 



.914 

-.0220 

-.0220 

-.0127 

-.0687 

-.0100 

-.0020 
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TABLE 12 


Plateau Values of Suction Ratio and Slopes of 
Suction-Ratio Segments - NACA 65A012 and 65,412 


CONFIGURATION 
A A 

h 

R 1 

d « 2 

F 

W 

AIRFOIL 




70 

45 

65A012 

Min. 

. 880 

0 

75 

45 

65A012 

1 

. 850 

-.0135 

80 

45 

6 5 A0 12 

i 

. 841 

-.0775 

70 

45 

65A012 

Std. 

.905 

0 

75 

45 

6 5 A0 12 

1 

.900 

-.0033 

80 

45 

65A012 

J 

. 860 

-.0254 

70 

45 

65A012 

Max. 

. 920 

0 

75 

45 

65A012 

T 

. 890 

-.0042 

80 

45 

6 5 A0 12 

1 

.880 

-.0222 

70 

45 

65]412 

Min. 

.930 

-.0047 

75 

45 

65 1 412 

J 

.900 

-.0080 

80 

45 

65 x 412 

f 

.925 

-.0156 

70 

45 

65 1 412 

Std. 

.930 

-.0047 

75 

45 

65i412 


.910 

-.0074 

80 

45 

65 x 412 

♦ 

.895 

-.0174 

70 

45 

65j412 

Max. 

.910 

0 

75 

45 

65,412 

T 

. 890 

-.0085 

80 

45 

65 1 412 

l 

.895 

-.0174 

45 

45 

00SLE 

Std. 

.968 

-.1003 

75 

45 

00SLE 

| 

.870 

-.1170 

80 

45 

00SLE 

1 

.890 

-.0988 


dR 

da 3 

dR 

dtt 4 

dR 

da 5 

dR 

da 6 

dR 

da 7 

-.0575 

-.0967 

-.0392 

-.0765 

-.0495 

-.0123 

-.0325 

-.0022 

-.0004 

-.0100 

+.0004 

-.0050 

-.0205 

0.0160 

-.0138 

-.0352 

-.0452 

-.0492 

-.0577 

-.0613 

-.0492 

-.0167 

-.00,12 

-.0024 

-.0167 

-.0012 

-.00054 

-.0058 
-.0180 
- . 0054 

-.0202 

-.0243 

-.0222 

-.0545 

-.0635 

-.0677 

-.0290 

-.0575 

-.0022 

-.0240 

-.0010 

-.0062 

-.0240 

-.0170 

-.0062 

-.0322 

-.0965 

-.0745 

-.1110 

-.0376 

-.0123 

-.0020 

-.0070 

-.0080 

-.0255 

-.0184 

-.0032 

-.0124 

-.0184 

-.0120 

-.0027 

-.0618 

-.0174 

-.0372 

-.0537 

-.0508 

-.0527 

-.0218 

-.0270 

-.0098 

-.0153 

-.0098 

-.0168 

-.0153 

-.0054 

-.0050 

-.0203 

-.0174 

-.0282 

-.0328 

-.0572 

-.0828 

-.0638 

-.0081 

-.0220 

-.0178 

-.0036 

-.0220 

-.0202 

-.0036 

-.0301 
-.0208 
- .0330 

-.0267 

-.0050 

-.0092 

-.0123 

-.0084 

-.0048 

-.0060 

-.0031 

-.0027 

-.0031 

-.0048 

-.0027 
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TABLE 13 PITCHING-MOMENT INTERCEPTS AND SLOPES OF PITCHING MOMENT-SEGMENTS - 

NACA OOXX AIRFOILS, RN^.^ 


CONFIGURATION 

C 

™a«o 

\,x 

C 

n 

a 2 

C 

m 

a 3 

C 

4 

C 

C 

C 

A 

F 

A 

w 

AIRFOIL 









45 

45 

0004 

-.008 

-.0019 

-.0030 

-.0013 

-.0042 

-.0010 

-.0037 

-.0010 



0008 

-.005 

-.0021 

-.0036 

-.0043 

-.0044 

+.0012 

-.0005 

-.0020 



0012 

-.003 

-.0014 

-.0030 

-.0030 

-.0047 

-.0028 

+.0002 

+.0006 



0015 

+.001 

-.0013 

-.0025 

-.0040 

-.0056 

+.0003 

-.0023 

-.0016 

60 

45 

0004 

-.007 

-.0034 

-.0050 

-.0026 

-.0035 

-.0010 

-.0017 

-.0024 

65 

45 

0004 

-.0065 

-.0048 

-.0058 

-.0038 

+.0001 

+.0001 

-.0021 

-.0021 



0008 

-.005 

-.0043 

-.0062 

-.0057 

-.0045 

-.0013 

+.0012 

-.0012 



0012 

-.002 

-.0039 

-.0059 

-.0082 

-.0054 

+.0020 

-.0002 

-.0008 



0015 

+.001 

-.0030 

-.0052 

-.0052 

-.0072 

+.0016 

-.0001 

-.0001 

70 

45 

0004 

-.003 

-.0056 

-.0073 

-.0056 

-.0056 

+.0037 

-.0029 

+.0040 



0008 

-.002 

-.0058 

-.0082 

-.0067 

-.0065 

+.0039 

-.0009 

+.0010 



0012 

-.002 

-.0050 

-.0068 

-.0078 

-.0019 

+.0017 

-.0010 

+.0010 



0015 

+.002 

-.0047 

-.0062 

-.0088 

-.0088 

+.0010 

+.0010 

+.0010 

75 

45 

0004 

-.006 

-.0070 

-.0084 

-.0076 

-.0118 

+.0062 

+.0007 

+.0007 



0008 

+.002 

-.0073 

-.0092 

-.0092 

-.0078 

-.0132 

+.0062 

-.0010 



0012 

+.005 

-.0065 

-.0084 

-.0084 

-.0106 

-.0043 

-.0108 

+.0076 



0015 

+.005 

-.0062 

-.0069 

-.0082 

-.0106 

-.0029 

-.0065 

+.0090 

80 

45 

0004 

+.001 

-.0087 

-.0098 

-.0108 

-.0109 

-.0126 

+.0106 

+.0106 



0008 

-.002 

-.0088 

-.0103 

-.0103 

-.0100 

-.0136 

-.0092 

+.0228 



0012 

+.001 

-.0083 

-.0108 

-.0108 

-.0078 

-.0132 

-.0096 

+.0137 



0015 

+.004 

-.0077 

-.0103 

-.0103 

-.0105 

-.0068 

-.0114 

-.0051 
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TABLE 14 PITCHING-MOMENT INTERCEPTS AND SLOPES OF PITCHING-MOMENT SEGMENTS - 

NACA OOXX AIRFOILS, RN___ 

STD 


CONFIGURATION 
^F AIRFOIL 

C 

m 

a«o 

C 

ra 

°fa.i 

C 

C 

in 

«3 

C 

C 

C 

C 

TD 

<*7 

45 

45 

0004 

-.008 

-.0019 

-.0030 

-.0013 

-.0042 

-.0010 

-.0037 

-.0010 



0008 

-.005 

-.0021 

-.0036 

-.0043 

-.0044 

+.0008 

+.0019 

-.0002 



0012 

-.003 

-.0014 

-.0033 

-.0033 

-.0054 

-.0008 

+.0007 

+.0004 



0015 

-.001 

-.0013 

-.0031 

-.0043 

-.0058 

+.0018 

+.0002 

0 

60 

45 

0004 

-.007 

-.0034 

-.0050 

-.0034 

-.0017 

-.0034 

-.0014 

-.0029 

65 

45 

0004 

-.0065 

-.0048 

-.0058 

-.0043 

+.0005 

-.0022 

-.0004 

-.0033 



0008 

-.005 

-.0043 

-.0062 

-.0077 

+.0010 

-.0009 

-.0009 

+.0005 



0012 

-.002 

-.0039 

-.0059 

-.0059 

-.0078 

0 

+.0028 

0 



0015 

+ .001 

-.0033 

-.0052 

-.0064 

-.0064 

-.0060 

+.0049 

0 

70 

45 

0004 

-.003 

-.0056 

-.0073 

-.0056 

-.0090 

+.0058 

-.0027 

+.0005 



0008 

-.002 

-.0058 

-.0082 

-.0070 

-.0054 

-.0077 

0 

0 



0012 

-.002 

-.0050 

-.0068 

-.0088 

-.0022 

-.0043 

+.0080 

0 



0015 

+.002 

-.0047 

-.0064 

-.0087 

-.0087 

-.0019 

-.0015 

+.0088 

75 

45 

0004 

-.006 

-.0070 

-.0084 

-.0077 

-.0114 

+.0062 

+.0007 

+.0007 



0008 

.002 

-.0073 

-.0077 

-.0095 

-.0063 

-.0128 

+.0120 

+.0012 



0012 

.005 

-.0065 

-.0084 

-.0084 

-.0106 

-.0010 

-.0106 

+.0126 



0015 

.005 

-.0062 

-.0069 

-.0082 

-.0106 

-.0029 

-.0065 

+.0063 

80 

45 

0004 

.001 

-.0087 

-.0098 

-.0108 

-.0108 

-.0126 

+.0106 

+.0106 



0008 

-.002 

-.0088 

-.0103 

-.0103 

-.0100 

-.0141 

-.0117 

-.0308 



0012 

.001 

-.0083 

-.0103 

-.0103 

-.0117 

-.0082 

-.0123 

+.0185 



0015 

.004 

-.0077 

-.0103 

-.0103 

-.0110 

-.0075 

-.0102 

-.0132 
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TABLE 15 PITCHING- MOMENT INTERCEPTS AND SLOPES FOR PITCHING-MOMENT SEGMENTS - 

NACA OOXX AIRFOILS, RN 

MAX 


CONFIGURATION 

C 

C 

C 

C 

C 

C 

C 

C 


A 

AIRFOIL 

P 3 

H 

O 

m C¥ 

0,1 

m 

<X 2 


m a 

4 


m 

«6 

m 

°7 

F 

w 








45 

45 

0004 

-.008 

-.0019 

-.0030 

-.0013 

-.0042 

-.0010 

-.0037 

-.0010 



0008 

-.005 

-.0021 

-.0036 

-.0049 

-.0049 

-.0038 

+.0006 

+.0015 



0012 

-.003 

-.0014 

-.0030 

-.0053 

-.0045 

+ .0019 

+.0008 

+.0002 



0015 

-.001 

-.0013 

-.0031 

-.0043 

-.0053 

+.0065 

+.0010 

-.0008 

60 

45 

0004 

-.007 

-.0034 

-.0050 

-.0034 

-.0017 

-.0034 

-.0014 

-.0029 

65 

45 

0004 

-.0065 

-.0048 

-.0058 

-.0043 

+.0005 

-.0034 

-.0002 

-.0033 



0008 

-.005 

-.0043 

-.0062 

-.0041 

-.0023 

-.0002 

-.0017 

0 



0012 

-.002 

-.0039 

-.0059 

-.0067 

-.0064 

-.0059 

+.0008 

+.0018 



0015 

+.001 

-.0036 

-.0049 

-.0049 

-.0060 

-.0036 

-.0011 

+.0018 

70 

45 

0004 

-.003 

-.0056 

-.0073 

-.0057 

-.0104 

+.0073 

-.0027 

+.0005 



0008 

-.002 

-.0058 

-.0053 

-.0053 

-.0051 

+.0030 

+.0030 

-.0005 



0012 

-.002 

-.0050 

-.0068 

-.0099 

-.0032 

-.0048 

+.0045 

0 



0015 

+.002 

-.0047 

-.0064 

-.0092 

-.0158 

-.0022 

+.0050 

+.0018 

75 

45 

0004 

-.006 

-.0070 

-.0084 

-.0077 

-.0114 

+.0062 

-.0070 

-.0070 



0008 

+.002 

-.0073 

-.0077 

-.0086 

-.0078 

-.0133 

+.0113 

+.0005 



0012 

+.005 

-.0065 

-.0084 

-.0093 

-.0054 

-.0054 

-.0100 

+.0077 



0015 

+.005 

-.0062 

-.0073 

-.0097 

-.0054 

-.0065 

-.0065 

+.0072 

80 

45 

0004 

+.001 

-.0087 

-.0098 

-.0108 

-.0108 

-.0126 

+.0106 

+.0106 



0008 

-.002 

-.0088 

-.0103 

-.0103 

-.0100 

-.0136 

-.0120 

+.0220 



0012 

+.001 

-.0083 

-.0105 

-.0102 

-.0106 

-.0132 

-.0098 

-.0040 



0015 

+.004 

-.0077 

-.0103 

-.0103 

-.0113 

-.0045 

-.0108 

-.0078 



TABLE 16 PITCHING- MOMENT INTERCEPTS AND SLOPES OF PITCH- 
ING-MOMENT SEGMENTS - NACA 65A012 AND 65j412 


CONFIGURATION 
^F A W AIRFOIL 

1 

Si 

C 

“a-o 

C 

■*0.1 

C 

a 2 

c 

"a 3 

C 

“<*4 

C 

m a 5 

C 

m <* 6 

c 

70 

45 

65A012 

MIN. 

+.005 

-.0049 

-.0063 

-.0048 

+.0025 

+.0062 

+.0005 

+.0023 

75 

45 

65A012 



+.003 

-.0058 

-.0070 

-.0087 

-.0059 

-.0072 

-.0037 

+.0115 

80 

45 

65A012 

1 


+.004 

-.0079 

-.0100 

-.0080 

-.0108 

-.0115 

-.0042 

+.0153 

70 

45 

65A012 

S'] 

fD. 

+.005 

-.0048 

-.0063 

-.0072 

-.0056 

-.0018 

+.0065 

0 

75 

45 

65 AO 12 



+.006 

-.0058 

-.0076 

-.0080 

-.0025 

-.0076 

-.0025 

+.0123 

80 

45 

65A012 



+.004 

-.0079 

-.0100 

-.0100 

-.0077 

-.0128 

-.0055 

+.0156 

70 

45 

65A012 

MAX 

+.005 

-.0048 

-.0063 

-.0060 

-.0093 

-.0039 

+.0069 

+.0069 

75 

45 

65 AO 12 



+.006 

-.0058 

-.0072 

-.0070 

-.0088 

0 

-.0074 

+.0175 

80 

45 

65A012 

1 


+.004 

-.0078 

o 

0 

o 

1 

-.0100 

-.0067 

-.0134 

-.0067 

+.0150 

70 

45 

65j412 

MIN. 

-.054 

-.0057 

-.0074 

-.0075 

-.0036 

-.0036 

+.0042 

+.0042 

75 

45 

65^412 



-.050 

-.0068 

-.0072 

-.0063 

-.0051 

-.0034 

.0145 

+.0081 

80 

45 

65j412 



-.045 

-.0087 

-.0096 

-.0102 

-.0100 

-.0110 

-.0066 

+.0123 

70 

45 

65 j4 12 

SI 

fD. 

-.053 

-.0053 

-.0072 

-.0072 

-.00835 

-.0011 

-.0057 

0 

75 

45 

65^12 



-.048 

-.0066 

-.0076 

-.0092 

-.0064 

0 

+.0156 

+.0077 

80 

45 

65^12 



-.045 

-.0087 

-.0096 

-.0096 

-.0110 

-.0082 

-.0115 

+.0057 

70 

45 

65j412 

MAX. 

-.053 

-.0053 

-.0072 

-.0074 

-.0056 

-.0056 

+.0048 

+.0048 

75 

45 

65j4l2 



-.048 

-.0066 

-.0076 

-.0066 

-.0092 

-.0058 

-.0004 

+.0162 

80 

45 

65j412 

\ 


-.045 

-.0087 

-.0096 

-.0096 

-.0106 

-.0085 

-.0079 

+.0095 

45 

45 

OOSLE 

STD. 

+.003 

-.0010 

-.0025 

-.0012 

-.0048 

+.0003 

+.0016 

-.0032 

75 

45 

OOSLE 



+.004 

-.0068 

-.0080 

-.0088 

+.0113 

+.0082 

+.0005 

-.0028 

80 

45 

OOSLE 



-.004 

-.0080 

-.0105 

-.0097 

-.0109 

-.0136 

-.0008 

+.0124 
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Figure 176. Effect of Airfoil-Thickness Ratio for Irregular 
Planform From LTPT Facility 
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form From LTPT Facility 










Figure 181. Effects of Airfoil Thickness Distribution and Camber 
for 75-45 Irregular Planform 






70-45 Irregular Planform 







Figure 184. Effect of Reynolds Number for NACA 65A012 Airfoil 






Figure 186, Effect of Fillet Sweep on Lift-Parameter Increment 
for Thin Airfoil With Sharp Leading Edge 
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Correlation for NACA 0004 Airfoils 













































































Lift Correlation for NACA 65A012 and 0012 Airfoils 
at Standard Reynolds Number 






Figure 203. Incremental Values of Lift 
for NACA 65A012 Airfoils a 
Reynolds Number 



0 

for 

32° 

0 

for 

30° 

0 

for 

28° 

0 

for 

26° 

0 

for 

24° 

0 

for 

22° 

0 

for 

20° 

0 

for 

O 
00 
1 — 1 

0 

for 

16° 

0 

for 

14° 

0 

for 

12° 

0 

for 

o 

o 

H 

0 

for 

8° 

0 

for 

6° 



Parameter 

Minimum 




0 for 20 
0 for 18 

0 for 16 
0 for 14 

0 for 12 


Figure 204. Incremental Values of Lift Parameter 
for NACA 65A012 Airfoils at Maximum 
Reynolds Number 
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Figure 208. Lift Correlations for Thin Airfoils 








mini 



32 ' 


E 


Figure 209. Lift Correlation of Wing III Planforms with NACA 
0008 Airfoils at Standard Reynolds Number from 
Test ARC 12-086 
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COMPARISONS WITH TEST DATA 


The only evaluation that could be made of the prediction 
methods described in the previous sections was to compare- predic- 
tions with test data from which the methods had been derived. In 
the comparisons which are presented in this section, the predic- 
tions were obtained using a computer program written for a CYBER 
172 computing system. The computer procedure utilizes simple 
linear interpolation between tabular data representations of many 
of the charts and graphs presented in this report. Therefore, 
one of the objectives of the comparisons of predictions with test 
data was to evaluate the accuracy of the computer procedure. A 
second objective was to evaluate what improvements had been 
achieved over the capabilities to predict the lift, drag, and pitch 
ing moment that had existed prior to this work. 

The first set of comparisons used the "basic" prediction 
method applied to the 9 configurations which were presented 
earlier in Figures 35 through 43 of this report as representative 
of the evaluation of the aeromodule prediction methods. Table 17 
identifies the specific figures which present these comparisons. 


Table 17 Key to Comparisons Made using the 
"Basic" SHIPS Prediction Method 


CONFIGURATION 


FIGURE 

NUMBER 


a f 

A w 

Airfoil 

Cl vs Of 

Cl vs Cd 

Cjh vs a 

Cm vs C L 

25 

25 

0008 

211a 

211b 

211c 

211d 

45 

45 

0008 

212a 

212b 

212c 

212d 

60 

60 

0008 

- 213a 

213b 

213c 

213d 

80 

25 

0008 

214a 

214b 

214c 

214d 

1 80 

45 

0008 

215a 

215b 

215c 

215d 

80 

60 

0008 

216a 

216b 

216c 

216d 

35 

25 

0008 

217a 

217b 

217c 

217d 

55 

45 

0008 

218a 

218b 

218c 

218d 

65 

45 

0008 

219a 

219b 

“219c 

219d 
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Figures 211, 212 and 213 present comparisons of predictions 
with test data from test ARC 12-086 for the basic Wing configura- 
tions Wi, Win, and Wv» respectively. The comparable comparisons 
with Aeromodule methods were shown in Figures 35, 36, and 37. For 
each of these configurations, the basic SHIPS prediction method 
produces much better agreement with test data than the Aeromodule 
methods for lift and pitching- moment variations with angle of at- 
tack and for drag variations with lift coefficient. Also shown 
for the SHIPS method are comparisons of predicted and test varia- 
tions of pitching moment with lift coefficient. Despite the fact 
that some discrepancies occur in the pitching -moment curves, the 
general character of the variations is well represented. 


Figures 214, 215 and 216 present similar comparisons for the 
irregular planforms consisting of Wings I, III, and V with 80 de- 
gree fillets. Corresponding comparisons with Aeromodule predic- 
tions were shown in Figures 38, 39 and 40. For these configura- 
tions the agreement of the SHIPS predictions with test data is 
very good for all of the characteristics and again is better 
than the Aeromodule methods. 

Figures 217, 218, and 219 present comparisons for irregular 
planforms consisting of Wings I, IX and III with fillet sweeps 
of 35, 55, and 65 degrees, respectively. Corresponding comparisons 
of Aeromodule predictions with test data were presented in Fig- 
ures 41, 42, and 43. Again, the SHIPS prediction method produces 
good agreement with the test data and much better agreement than 
the Aeromodule method for all of the characteristics. 

In general, the "basic" SHIPS prediction method produces good 
agreement with test data and represents a significant improvement 
over the Aeromodule method from which some of the components of 
the SHIPS method were adapted. There is one disturbing trend 
that should be noted in the SHIPS predictions. At the high an- 
gles of attack, there is a tendency to predict the lift coefficient 
slightly lower -or slightly higher than the test data. These 
slight discrepancies are reflected in the pitching- moment varia- 
tions with angle of attack and sometimes produce an exaggerated 
or even erroneous change in slope. This situation occurs because 
the stall progression function used in the lift prediction' 
does not completely account for all of the variatiQns 
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which occur due to outboard panel sweep variations. As a conse- 
quence, the user should pay attention to the pitch-up/pitch-down 
design guides when evaluating a specific configuration. 

Ten additional comparisons of predictions with test data are 
presented, some of which use only the "basic” method, some of 
which use both the "basic" and the "revised" methods, and some of 
which use only the "revised" method. Table 18 presents a key to 
these additional comparisons with test data. 

Figure 220 presents comparisons of predictions for an ir- 
regular planform having a fillet sweep of 75 degrees, an outboard 
panel sweep of 40 degrees, and an NACA 0008 airfoil. The outboard 
panel sweep is intermediate between two of the test configura- 
tions, and the basic wing planform fits this SHIPS planform family 
in terms of half-chord sweep, taper ratio, and aspect ratio. The 
test data show that, except for the variation of suction ratio 
with angle of attack, the test data bracket the predictions. 

The discrepancy for the suet ion -ratio curve represents one of the 
compromises that had to be made during prediction method develop- 
ment to achieve smooth variations of the angle-of-attack bound- 
aries and slopes in the quasi-linear representation. 

Figure 221 presents comparisons of prediction with test for 
the irregular planform configuration 60 25.0008 using the basic 
SHIPS prediction methods. The comparison for lift, suction ratio, 
and drag are very good. The variations of pitching moment with 
angle of attack and lift coefficient both show some discrepan- 
cies. The discrepancies for the irregular planforms with rela- 
tively low values of fillet sweep and outboard panel sweep 
occur because the linear segment representations of 
the pitching -moment variation with angle of attack requires an 
additional segment which is not reflected in the suction- ratio 
variation with angle of attack. The "basic" SHIPS method utilized 
the angle-of-attack boundaries determined from the suction ratio 
analysis, and provisions were made for use of fictitious values to 
accommodate the pitching-moment curves only at low and very high 
angles of attack. 

Figure 222 presents comparisons of predictions made with 
both the "basic" and "revised" SHIPS methods and test data from 
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Table 18 



229a 229b 229c 229d 229e 


















two different facilities but obtained at the same Reynolds 
number (26.25 million per meter). Ihe predictions made using 
the "basic" method slightly undershoot the lift curve at high 
angles of attack which is reflected in the drag polar and 
pitching -moment variation with lift coefficient for the 
irregular planform 65 45.0008. The "revised" prediction method 
very accurately predicts all three characteristics. The 
differences between the sets of test data are representative 
of the magnitudes which occurred for all of the configurations 
which were tested in both facilities. This is a representation 
of the uncertainties in wind-tunnel test data historically 
noted for configurations tested in several facilities. 

Figures 223, 224 and 225 present data comparisons with 
predictions for Basic wing III with NACA 0004, NACA 0008 and 
NACA 0012 airfoil sections, respectively. The "revised" 
prediction methods produce accurate representations of the 
test data for all three airfoil thickness ratios with the 
exception of the drag at low coefficients for the NACA 0004 
airfoil configuration. The measured minimum drag for that 
configuration appeared to be unusually low at the test Mach 
number of 0.20, probably due to balance accuracy. 

Figure 226 shows the comparisons of predictions using the 
"revised" prediction method with test data for the irregular 
planform configuration 75 45.0004. For this case, the method 
predicts the test data quite accurately even though the pre- 
diction was inadvertently made for Mach 0.30 and the data were 
obtained at M = 0.20. Note that the test minimum drag value 
is the same for this configuration as it was for the 45 45.0004 
configuration which is further evidence that the test value for 
the 45 45.0004 configuration is unusually low. Of particular 
importance is the fact that the pitch-up at high angle of 
attack is accurately predicted. 

Figure 227 presents comparisons of predictions with test 
data for irrjegular planform 75 45.0012 for the standard and 
maximum Reynolds number conditions. The revised prediction 
methods produce good agreement for both Reynolds number con- 
ditions and, in particular, reflect the increase in angle of 
attack for maximum lift coefficient and subsequent pitch-up. 
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The next comparison of predictions with test data is for the 
irregular planform consisting of Basic Wing III with 75 degree 
sweep fillet and having the cambered NACA 65,412 airfoil and is 
presented in Figure 228. In this case the effect of camber on 
the lift, drag, and pitching moment curves are well predicted with the ex- 
ception of a slight undershoot on the lift curve at the highest 
angles of attack which is reflected in the drag polar and pitch- 
ing moment curves. Again the shape of the predicted pitching 
moment curve reflects the pitch-up characteristic which occurs 
in the test data. 

The last comparisons between predictions and test data 
presented in Figure 229 use both the ’'basic" and "revised" pre- 
diction methods and are for the irregular planform consisting 
of Wing III with the 80 degree fillet having the NACA 0008 airfoil. 
These comparisons show that the "revised" method more accurately 
predicts the lift curve than does the "basic" method. This 
improvement is also shown to correct an erroneous occurrence in 
the pitching -moment variation with lift coefficient in which the 
"basic" method shows a pitch- down whereas the "revised" method 
accurately predicts a pitch-up. In retrospect it appears that 
the simplified representation of the lift curve in the "basic" 
prediction method using a single curve for the stall progression 
function should be modified to better account for the effects 
of outer panel planform effects. The "revised" method is limited 
in the scope of planforms to which it can be applied, but pro- 
duces very accurate results and is more applicable when predict- 
ing full-scale conditions of airfoils having higher thickness 
ratios. 

The fact that the test results obtained in different facili- 
ties are somewhat different reflects the fact that even at low 
speed there are still uncertainties inherent in the available 
test techniques. 



(a) Lift Curve 

211. Comparisons Between Basic SHIPS Prediction and Test 
for Basic Wing I 





(b) Drag Polar 
Figure 211 . Continued 
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(c) Pitching-Moment Variation with Angle of Attack 
Figure 211. Continued 
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(d) Pitching -Moment Variation with Lift 
Figure 211. Concluded. 








2X2. Comparisons Between Basic SHIPS Prediction and Test 
for Basic Wing III 
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212. Continued 





(c) Pitching -“Moment Variation with Angle of Attack 
Figure 212. Continued 
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(d) Pitching-Moment Variation with Lift 
gure 212. Concluded. 
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Figure 213. Continued 





(c) Pitching-Moment Variation with Angle of Attack 
Figure 213. Continued. 
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(d) Pitching-Moment Variation with Lift 
Figure 213. Concluded. 
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(a) Lift Curve 

Figure 214. Comparisons Between Basic SHIPS Prediction and Test 
Wing I with 80° Fillet 







Figure 2X4. Continued 











(d) Pitching-Moment Variation with Lift 
Figure 214. Concluded 








(a) Lift Curve 

Figure 215. Comparisons Between Basic SHIPS Prediction and Test 
Wing III with 80° Fillet 






(b) Drag Polar 
Figure 215. Continued 







Test 086-83 
X6I Basic 
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(c) Pitching-Moment Variation with Aigle of Attack 
Figure 215. Continued 








(d) Pitching-Moment Variation with Lift 
Figure 215. Concluded. 










Figure 216. Continued 
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(c) Pitching-Moment Variation with Angle of Attack 
Figure 216. Continued 
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(d) Pitching -Moment Variation with Lift 
ire 216. Concluded. 






(a) Lift Curve 

Figure 217. Comparisons of Basic SHIPS Prediction and Test 
for Basic Wing I with 35° Fillet 






Figure 217. Continued 
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(c) Pitching -Moment Variation with Angle of Attack 





(d) Pitching-Moment Variation with Lift 
re 217. Concluded. 







(a) Lift Curve 

218. Comparisons of Basic SHIPS Prediction and Test 
for Wing III with 55° Fillet 







(b) Drag Polar 
Figure 218. Continued 





(c) Pitching-Moment Variation with Angle of Attack 
Figure 218. Continued 






(d) Pitching-Moment Variation with Lift 
re 218. Concluded. 
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(a) Lift Curve 

Figure 219. Comparisons of Basic SHIPS Prediction and Test 
for Wing V with 65° Fillet 
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(b) Drag Polar 
Figure 219. Continued 





(c) Pitching -Moment Variation with Angle of Attack 






(d) Pitching-Moment Variation with Lift 
Figure 219. Concluded 
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(c) Drag Polars 
Figure 220. Continued 













(d) Pitching-Moment Variation with Angle of Attack 
Figure 220. Continued 

















(e) Pitching- Moment Variations with Lift 
Figure 220. Concluded. 













(a) Lift Curve 

Figure 221. Comparisons of Basic SHIPS Prediction and Test 
for Wing I with 60° Fillet 


















Figure 221. Continued 







(c) Drag Polar 
Figure 221. Continued 
























(d) Pitching-Moment Variation with Angle of Attack 
























(e) Pitching-Moment Variation with Lift 
e 221. Concluded 
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Figure 222. Comparisons of Basic and Revised SHIPS Predictions 
and Test for Wing III with 65° Fillet for Two 
Different Test Facilities 










Figure 222. Continued 



















(d) Pitching -Moment Variation with Lift 
Figure 222. Concluded 















(a) Lift Curve 

Figure 223. Comparisons of Revised SHIPS Prediction and Test 
for Wing III with NACA 0004 Airfoil 
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(b) Suction-Ratio Variation with Angle of Attack 
Figure 223. Continued 













Figure 223. Continued 




























(e) Pitching Moment Variation with Lift 
Figure 223. Concluded 







(a) Lift Curve 

Figure 224. Comparisons of Basic and Revised SHIPS Predictions 

and Test for Basic Wing III with NACA 0008 Airfoils for 
Two Different Facilities 















(c) Pitching-Moment Variation with Angle of Attack 
Figure 224. Continued 






(d) Pitching-Moment Variations with Lift 
Figure 224. Concluded 








(a) Lift Curve 

Figure 225. Comparisons of Revised SHIPS Prediction and Test 
for Wing III with NACA 0012 Airfoils 








'Hi Airfoil 
45 0012 



(b) Drag Polar 
Figure 225. Continued 











(c) Pitching-Moment Variation with Angle of Attack 
Figure 225. Continued 
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(d) Pitching -Moment Variation with Lift 
Figure 225. Concluded 









(a) Lift Curve 

Figure 226. Comparisons of Revised SHIPS Prediction and Test 

for Wing III with 75° Fillet and NACA 0004 Airfoils 






) Test 266-60 M 
- X6I Revised M, 



(b) Drag Polar 
Figure 226. Continued 













(c) Pitching -Moment Variation with Angle of Attack 
Figure 226. Continued 
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O Test 266 
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(d) Pitching Moment Variation with lift 
Figure 226. Concluded 








Airfoil 
45 0012 



(a) Lift Curve 

Figure 227. Comparisons of Revised SHIPS Predictions and Test 

for Wing III with 75° Fillet and NACA 0012 Airfoils 
Show Effects of Reynolds Number 



















'W Airfoil 
45 0012 



(b) Drag Polar 
Figure 227. Continued 











(c) Pitching-Moment Variation with Angle of Attack 
Figure 227. Continued 
















F Airfoil 



(a) Lift Curve 

Figure 228. Comparisons of Revised SHIPS Prediction and Test 
for Wing III with 75° Fillet and Cambered NACA 
65-.412 Airfoil 













(b) Drag Polar 
Figure 228. Continued 
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(c) Pitching- Moment Variation with Angle of Attack 
Figure 228. Continued 





(d) Pitching-Moment Variation with Lift 
Figure 228. Concluded 
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Figure 229. Comparisons of Basic and REvised SHIPS Predictions 
and Test for Wing III with 80° Fillet and NACA 0008 
Airfoils for Two Different Facilities 






(b) Suction-Ratio Variations with Angle of Attack 
Figure 229. Continued 



















(c) Drag Polar 
Figure 229 . Continued 





























(e) Pitching-Moment Variations with Lift 
Figure 229. Concluded 


















CONCLUDING REMARKS 


The objectives of this investigation were to develop predic- 
tion methods for low speed lift, drag, and pitching moment of ir- 
regular planform wings and design guidelines appropriate for pre- 
liminary design of advanced aerospace vehicles. 

At the outset of the investigation it was hoped that quite 
simple empirical methods could be evolved that would be suitable 
for solution with nothing more than a hand calculator. As the 
investigation progressed, it became apparent that the flow fields 
produced by the component panels of irregular planform wings have 
complicated interactions that are affected by differences in Rey- 
nolds number in different ways depending on the fillet/wing com- 
bination and on the airfoil thickness. As a consequence the 
methods as finally devised are more complex than originally en- 
visioned and require a computer procedure for efficient solution. 

The methods provide a significant improvement in capability 
to accurately predict variations of lift, drag, and pitching- mo- 
ment with angle of attack for a wide range of irregular planforms, 
and for a more limited range of planforms can account for effects 
of airfoil -thickness ratio and Reynolds number. 

The investigation also showed that there are differences in 
test results obtained in different test facilities that introduce 
inherent uncertainties in the low speed aerodynamic characteris- 
tics of any given configuration. The magnitude of the uncertain- 
ties can be evaluated for the limited range of irregular plan- 
forms by making predictions using both the "Basic'’ and "Revised" 
prediction methods. 
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APPENDIX A 

REYNOLDS NUMBER EFFECT ON LIFT OF SHIPS WINGS 


During the investigation described in the main body of this 
report, a concerted effort was made to derive some simple method 
to obtain corrections to the lift curves for the effects of Rey- 
nolds number. No successful method was found that could ade- 
quately describe the effects as a general function of Reynolds 
number for the complete range of SHIPS planforms. 

In order to facilitate extrapolation of test results from 
a test performed at relatively low Reynolds numbers to flight 
conditions at low speed, data on the effects of Reynolds number 
on the lift correlation parameter are presented in this appendix 
for all configurations having airfoils with thickness ratios of 
0.08 or larger. There were no significant effects for the 
NACA 0004 or thin hexagonal airfoils. 

The data are presented in the form of the variation of the 
incremental value of the lift correlation parameter, ^/Cl Ai \ 

\^V 

with angle of attack for various values of unit Reynolds number. 
The baseline for the incremental value is the value of the lift 
correlation parameter obtained at the standard unit Reynolds num- 
ber of 26.25 million per meter at each angle of attack for each 
configuration. 

For the reader’s convenience the plots are sorted into two 
major categories: 

1. Corrections that were obtained from data for test 
ARC 12-086. These data should only be applied to 
predictions made using the "basic” prediction 
method (Figures 230, 231). 
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2. Corrections that wete obtained from Langley LTPT Tests 
255, 262 and 266. These data apply only to predic- 
tions made using the "revised" prediction method for 
the "standard" Reynolds number (Figures 232-235) . 

The appropriate values of Ct > A, and are listed in 
Table 19. 
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TABLE 19 NUMERICAL VALUES OF ELEMENTS IN LIFT CORRELATION PARAMETER 
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NOTE: T) b = .41157 for all wings 
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Figure 230. Incremental Values of Lift Correlation Parameter 

at Unit Reynolds Numbers of 13.13 and 19.67 Million 
Per Meter From Test ARC 12-086 for SHIPS Planforms 
Having Constant Values of Fillet Sweep 
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Figure 230. Continued 






Figure 230. Continued 






Figure 230. Continued 
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Figure 230. Continued 



Figure 230. Concluded 








421 


Figure 231. Incremental Values of Lift Correlation Parameter 

at Unit Reynolds Numbers of 13.13 and 19.67 Million 
Per Meter from Test ARC 12-086 Showing Effect of 
Thickness Ratio for Wing I With Various Fillet 
Sweeps 
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Figure 231. Continued 







Figure 231. Concluded 
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Per Meter from Various Langley LTPT Tests Showing 
Effect of Thickness Ratio for Wing III with Various 
Fillet Sweeps 



Figure 232. Continued 





Figure 232. Continued 







Figure 232. Continued 





gure 232. Concluded 
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at Unit Reynolds Numbers of 32.81, 39.37 and 45.93 
Million Per Meter from Various Langley LTPT Tests 
Showing Effects of Thickness Ratio for Wing III 
with Various Fillets 






Figure 233. Continued 
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Figure 233. Concluded 
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Per Meter from Langley LTPT Test 262 Showing Effects 
of Airfoil Section and Camber for Wing III with 
Various Fillet Sweeps 
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Figure 235. Incremental Values of Lift Correlation Parameter ai 
Unit Reynolds Numbers of 32.81, 39.37 and 45.93 
Million Per Meter from Langley LTPT Test 262 Show- 
ing Effects of Airfoil Section and Camber for Wing 
III with Various Fillet Sweeps 
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Figure 235. Concluded 
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